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ABSTRACT 

We present self-similar solutions that describe the gravitational collapse of rotating, isothermal, magnetic 
molecular-cloud cores. These solutions make it possible, for the first time, to study the formation of rotationally 
supported protostellar disks of the type detected around many young stellar objects in the context of a realistic sce- 
nario of star formation in magnetically supported, weakly ionized, molecular cloud cores. This work focuses on 
the evolution after a point mass first forms at the center and generalizes previous results by Contopoulos, Ciolek, 
& Konigl that did not include rotation. Our semianalytic scheme incorporates ambipolar diffusion and magnetic 
braking and allows us to examine the full range of expected behaviors and their dependence on the physical pa- 
rameters. We find that, for typical parameter values, the inflow first passes through an ambipolar-diffusion shock 
(at a radius r a ), where the magnetic flux decouples from the matter, and subsequently through a centrifugal shock 
(at r c ), inward of which a rotationally supported disk (of mass Mj) is established. By the time (~ 10 5 yr) that 
the central mass M c grows to ~ 1 M Q , r a > 10 3 AU, r c > 10 2 AU, and Mj /M c < 0. 1 . The derived disk properties 
are consistent with data on T Tauri systems, and our results imply that protostellar disks may well be Keplerian 
also during earlier phases of their evolution. We demonstrate that the disk is likely to drive centrifugal outflows 
that transport angular momentum and mass, and we show how the radially self-similar wind solution of Bland- 
ford & Payne can be naturally incorporated into the disk model. We further verify that gravitational torques and 
magnetorotational instability-induced turbulence typically do not play an important role in the angular momentum 
transport. For completeness, we also present solutions for the limiting cases of fast rotation (where the collapse 
results in a massive disk with such a large outer radius that it traps the ambipolar-diffusion front) and strong 
braking (where no disk is formed and the collapse resembles that of a nonrotating core at small radii), as well as 
solutions describing the rotational collapse of ideal-MHD and of nonmagnetic model cores. 

Subject headings: accretion, accretion disks — diffusion — ISM: clouds — ISM: magnetic fields — MHD — 
shock waves — stars: formation 



1. INTRODUCTION 

Low-mass star formation is believed to occur predominantly 
in molecular cloud cores. The commonly accepted scenario is 
that dense molecular clouds are supported against self-gravity 
by magnetic and thermal stresses. Since the gas is only weakly 
ionized, the clouds are not in a strict steady state: the prevalent 
neutral molecules, pulled in by the force of gravity, gradually 
drift inward through the magnetic field (attached to the ionized 
component and anchored in the cloud envelope) in a process 
known as ambipolar diffusion. When the mass-to-flux ratio 
in the cloud core exceeds a critical value, the core undergoes 
gravitational collapse to form a central protostar. The core evo- 
lution after the collapse is initiated can in general be divided 
into two phases: a rapid initial dynamical contraction, which 
proceeds nonhomologously and eventually results in the for- 
mation of a central point mass, followed by an accretion phase, 
during which the central mass gradually increases and its grav- 
itational field dominates a progressively larger region around 
the origin. The evolution of star-forming molecular clouds has 
been studied mostly through numerical simulations (see, e.g., 
Mouschovias & Ciolek 1999 for a review). Until recently, these 
calculations were terminated before the end of the dynamical- 
collapse phase, basically at the point where radiative trapping 
near the center started to invalidate the isothermality assump- 
tion adopted in these computations. However, Ciolek & Konigl 
(1998, hereafter CK) circumvented this problem by treating the 
comparatively small region of radiative trapping as a central 



sink cell; this enabled them to model the collapse of nonrotating 
magnetic cloud cores through point-mass formation (PMF) and 
into the protostellar accretion phase. Concurrently, Contopou- 
los, Ciolek, & Konigl (1998, hereafter CCK) obtained semi- 
analytic self-similar solutions that explicitly accounted for the 
effects of ambipolar diffusion and successfully reproduced the 
main qualitative features of the numerical simulations of CK. 

Our goal in this paper is to extend the previous work on 
the post-PMF evolution of collapsing magnetic cloud cores 
by including the effects of rotation. Molecular line ob- 
servations (e.g., Goodman et al. 1993; Kane & Clemens 
1997) have established that a majority of dense (> 10 4 cm -3 ) 
cloud cores show evidence of rotation, with angular velocities 
~ 3xl0 _15 -10 _13 s _1 that tend to be uniform on scales of 
~ 0.1 pc, and with specific angular momenta in the range 
~ 4x 10 20 -3x 10 22 cm 2 s _1 . Although present, the rotation con- 
tributes only a small fraction (typically no more than a few per- 
cent) to the dynamical support of these cores, and self-gravity is 
mostly balanced by magnetic and thermal stresses. During their 
quasistatic contraction phase, the cores evidently lose angu- 
lar momentum by means of magnetic braking — the magnetic 
transfer of angular momentum to the ambient gas through tor- 
sional Alfven waves — and this process also tends to align their 
angular velocity and large-scale magnetic field vectors (e.g., 
Mouschovias & Ciolek 1999). Once dynamical collapse is ini- 
tiated and a core goes into a near-free-fall state, the specific 
angular momentum is expected to be approximately conserved, 
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resulting in a progressive increase in the centrifugal force that 
eventually halts the collapse and gives rise to a rotationally 
supported disk. These expectations are consistent with the 
results of molecular-line interferometric observations, which 
have probed contracting cloud cores on scales <~ 10 2 -10 3 AU 
(e.g., Myers, Evans, & Ohashi 2000; Mundy, Looney, & Welch 
2000; Wilner & Lay 2000). These observations have revealed 
that, on scales < 10 3 AU, the cores have a flattened, thick-disk 
morphology and a velocity field that is dominated by infall mo- 
tions (with only a modest rotational component). This morphol- 
ogy is consistent with numerical simulations of magnetically 
supported clouds, which have demonstrated that the gas rapidly 
contracts along the field lines and maintains force equilibrium 
along the field even during the collapse phase (e.g., Fiedler & 
Mouschovias 1993; Galli & Shu 1993; CK), including in cases 
where the clouds are initially elongated in the field direction 
(e.g., Nakamura, Hanawa, & Nakano 1995; Tomisaka 1996). 
The observations have also established that angular momentum 
is by and large conserved in the infalling gas (e.g., Ohashi et al. 
1997) and leads to the formation of rotationally supported disks 
on scales < 10 2 AU. 

The pre-PMF collapse of rotating, magnetized, weakly 
ionized cloud cores was previously studied by Basu & 
Mouschovias (1994, 1995a,b). They found that little angular 
momentum is lost during the dynamical-collapse phase but that 
the centrifugal force nevertheless does not become important 
prior to PMF. In the current work we focus on the post-PMF 
phase of the collapse, which enables us to model the formation 
of circumstellar disks in the context of magnetic cloud cores 
and to study their properties and their role in the buildup of 
the central protostar. Circumstellar disks have been detected in 
<~ 25-50% of pre-main-sequence stars in nearby dark clouds 
(e.g., Beckwith & Sargent 1993), and it is likely that most of 
the mass assembled in a typical low-mass young stellar object 
(YSO) is accreted through a disk (e.g., Calvet, Hartmann, & 
Strom 2000). The study of protostellar disks is important also 
because they are the incubators of planetary systems, so their 
properties are directly relevant to the process of planet forma- 
tion (e.g., Wuchterl, Guillot, & Lissauer 2000). The rotational 
properties of the collapsing core are evidently also of particular 
relevance to the question of multiple-star formation (e.g., Bo- 
denheimer et al. 2000). In the ensuing discussion we assume 
that the collapsing core does not give rise to more than one 
stellar object; we nevertheless briefly address the issue of core 
fragmentation in § 4. 

The studies of CK and CCK demonstrated that ambipolar 
diffusion, which is unimportant during the dynamical pre-PMF 
collapse phase, is "revitalized" after a point mass starts to grow 
in the center and leads to a decoupling of the magnetic flux 
from the inflowing gas. It was found that the "decoupling front" 
propagates outward in the form of a hydromagnetic shock (as 
originally suggested by Li & McKee 1996) and that this process 
can go a long way toward resolving the magnetic flux problem 
in star formation (the several-orders-of-magnitude discrepancy 
between the empirical upper limit on the magnetic flux of a 
protostar and the flux associated with the corresponding mass 
in the pre-collapse core). By incorporating rotation into the pic- 
ture, we can study how this issue is affected by the formation of 
a disk. Furthermore, as our model includes angular momentum 
transport (by the large-scale magnetic field), which is crucial 
for enabling mass to accumulate in the center and form a proto- 
star (see § 2.1), we are also in a position to address the angular 



momentum problem in star formation (analogous to the mag- 
netic flux problem, except that the quantity under consideration 
is the specific angular momentum). 

The emergence of a quasi-stationary, rotationally supported 
disk from a nearly freely falling, collapsing core inevitably 
involves a strong deceleration of the inflowing matter, which 
almost invariably takes place in a shock. This "centrifugal" 
shock is distinct from the ambipolar-diffusion shock mentioned 
above: it is typically hydrodynamic, rather than hydromag- 
netic, in nature, and it is generally located at a different dis- 
tance from the center. The appearance of a second shock in- 
creases the complexity of the problem, and it might seem that 
it could only be tackled through a multi-scale numerical simu- 
lation. However, we have found that it is possible to treat this 
problem semianalytically by generalizing the self-similar solu- 
tions presented in CCK to include the effects of rotation and 
magnetic braking. Although this treatment is somewhat less 
detailed than a numerical simulation, the semianalytic approach 
allows us to more readily explore the fairly extensive parame- 
ter space. This approach is supported by the strong evidence 
from numerical simulations that a multiscale core collapse nat- 
urally leads to a self-similar evolution. In particular, this be- 
havior has been found to describe the effectively magnetic flux- 
and angular momentum-conserving pre-PMF collapse phase of 
a rotating magnetic core seen in the aforementioned numerical 
simulations of Basu & Mouschovias (Basu 1997; Nakamura 
et al. 1999; see also Narita, Hayashi, & Miyama 1984 and Li 
& Shu 1997). As regards the post-PMF evolution, Saigo & 
Hanawa (1998) obtained similarity solutions for flattened cloud 
cores that approximate previous numerical simulations of the 
collapse of rotating clouds through the moment when a central 
point mass could be expected to form. The good correspon- 
dence between the CCK solutions and the nonrotating magnetic 
collapse simulations of CK further strengthens our confidence 
in the viability of this approach. 

In relating this work to previously published results, we note 
that neither the Saigo & Hanawa (1998) model nor the numeri- 
cal simulations with which they compare their solutions incor- 
porate an angular momentum transport mechanism, and there- 
fore they do not in practice give rise to a central point mass. 
(Nevertheless, as we discuss in § 3.1, their post-PMF solutions 
correspond to a limiting case of the more general model consid- 
ered in this paper.) It is also worth mentioning in this connec- 
tion that self-similar models of the gravitational collapse of ro- 
tating, disk-like mass distributions that incorporate viscous an- 
gular momentum transport have been constructed by Mineshige 
& Umemura (1997) and Tsuribe (1999). However, the for- 
mer model arbitrarily imposes an inner boundary condition that 
does not allow mass to accumulate at the center, whereas the 
latter model is restricted to the case in which the inflow speed is 
everywhere so low that no centrifugal shock forms. Neither of 
these models therefore provides a realistic framework for study- 
ing star formation in rotating molecular cloud cores. 

The plan of this paper is as follows. In § 2 we formulate the 
problem, discuss the adopted approximations, and describe the 
self-similar model. In § 3 we present illustrative solutions, or- 
dered by increasing level of complexity, which isolate different 
aspects of the full problem. We use the flexibility afforded by 
our semianalytic approach to examine the dependence of the 
results on the efficiency of ambipolar diffusion and magnetic 
braking and to elucidate the interplay between these two pro- 
cesses. We give a general discussion of these solutions and of 
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their astrophysical implications in § 4, where we also consider 
the production of centrifugally driven disk outflows. Our con- 
clusions are summarized in § 5. 

2. SELF-SIMILAR MODEL 
2. 1 . Formulation of the Problem 

Our basic model setup is the same as in CCK. We con- 
sider an isothermal, disklike core that contracts in the radial 
direction while maintaining a vertical hydrostatic equilibrium 
at all times. As noted in § 1, this behavior was found in nu- 
merical simulations of isothermal cores that are initially sup- 
ported against gravitational collapse by a large-scale, ordered 
magnetic field. Nonmagnetic simulations (e.g., Norman, Wil- 
son, & Barton 1980; Narita, Hayashi, & Miyama 1984; Mat- 
sumoto, Hanawa, & Nakamura 1997) have established that ro- 
tation could further contribute to the tendency of a core to flat- 
ten as it contracts so long as the specific angular momentum 
is approximately conserved. In restricting gas motions in our 
model to the plane of the flattened core, we preclude the possi- 
bility that matter can reach the center from above (or below) the 
disk, and this, in turn, implies that a point mass can form in the 
center only if there exists a mechanism for removing angular 
momentum from the inflowing gas. There have also been mod- 
els in which the complementary view has been adopted, viz., 
new mass is added to the protostar and its circumstellar disk 
only from above and below the disk plane (e.g., Terebey, Shu, 
& Cassen 1984; Yorke, Bodenheimer, & Laughlin 1993). These 
models tend to produce rotationally supported disks whose 
masses are comparable to that of the central star, even if an- 
gular momentum transport within the disk is taken into account 
(e.g., Cassen & Summers 1983; Lin & Pringle 1990). How- 
ever, most observed protostellar disks in nearby dark clouds 
have masses (estimated from dust emission or, when possible, 
from a direct measurement of the disk velocity field) that are 
significantly smaller (< 10%) than the protostellar mass (e.g., 
Beckwith & Sargent 1993). As we demonstrate in § 3, such 
low-mass systems can form naturally when mass flows into the 
central region primarily in the disk plane. Given that this is 
the expected situation for magnetically regulated core contrac- 
tion and that the disklike configurations predicted in this case 
are consistent with observations of collapsing cloud cores, we 
consider this approximation to be adequate for capturing the 
essence of protostellar-disk and star formation. 

The isothermality assumption was justified in the work of 
CK and CCK on the grounds that it generally applies so long 
as the density remains below <~ 10 10 cm~ 3 , which for a typical 
nonrotating collapse is the case on scales > 5 AU. This assump- 
tion is less tenable in the present work because, when a quasi- 
stationary, centrifugally supported disk forms, it has a signif- 
icantly higher column density at a given radius than a nearly 
freely falling configuration, with the result that radiative trap- 
ping already occurs on scales < 10 2 AU (e.g., Wardle & Konigl 
1993). Irradiation by the central protostar tends to mitigate this 
effect and establish vertical isothermality in the outer regions of 
the disk, but, in turn, it leads to a moderately strong (T oc r~' /2 ) 
variation of the disk temperature with radius (e.g., D' Alessio et 
al. 1998). However, since thermal stresses turn out not to play 
a major role in the core-collapse dynamics, the assumption of 
isothermality (which is required for self-similarity) is not likely 
to introduce significant inaccuracies into the results. 

Another requirement for self-similarity is that the ion density 

2 Our dimensionless ambipolar diffusion parameter r; = T„i (4ttGp) 1j/2 is smaller 



Pi be proportional to the square root of the neutral mass density 
(which is essentially equal to the total gas density p). In this 
case, the ratio rj of the neutral-ion momentum-exchange time 
scale r n i to the nominal self-gravity free-fall time (47rG/})~ 1/2 
is independent of density and can be treated as a constant. 2 This 
behavior is found at comparatively low densities, when molec- 
ular ionization by cosmic rays (at a rate £ = 10~ 17 £_i7S _1 ) is 
balanced by rapid dissociative recombination of the molecular 
ions and the dominant positive charges are metal ions (formed 
by a charge-exchange reaction between molecular ions and neu- 
tral metal atoms and destroyed by recombination on grain sur- 
faces). For grains of radius a = 0.1 pm, a temperature T = 10 K, 
and £_i7 w 1, this scaling applies roughly at neutral particle 
densities between <~ 10 4 cm~ 3 and < 10 7 cm -3 (e.g., Ciolek 
& Mouschovias 1998; Kamaya & Nishi 2000). These val- 
ues of T and n characterize the outer regions of collapsing 
cores, typically on scales > 10 3 AU (e.g., CK). In this regime, 

V ^ 0.2 C\i 2 - 

The numerical models of CK extended to scales < 10AU, 
where densities > 10 10 cm~ 3 were attained. In the case of a 
rotationally supported disk, even higher densities are expected 
on these scales. At such high densities, grains are typically the 
dominant carriers of both positive and negative charges, and 
their densities also scale as n 1 /2 (e.g., Nishi, Nakano, & Ume- 
bayashi 1991). Assuming, for the sake of illustration, that the 
charged grain distribution is dominated by a population of small 
(PAH-like) singly charged particles of radius a (see Neufeld & 
Hollenbach 1994), we can deduce the charged-grain number 
density n± from the ionization balance relation 

«»=*< 2 «> v *0" J ( i+ w)< <» 

where m± r = 2ira 3 p s /3 (with p s w 2.3 gem -3 ) is the reduced 
grain mass. Taking as fiducial values n = 10 12 «i2cm~ 3 and 
T = lO^K at a distance r w 10AU from the center of a ro- 
tationally supported protostellar disk, we can neglect the factor 
1 in comparison with the second term in the last parentheses 
on the right-hand side of equation (1) (which accounts for the 
electrostatic attraction between the recombining grains). We 
then get n± w 15.8(a/5 A) 1 / 4 7 , 2 1/4 £y i 2 7 nJ 2 ' 2 crn _3 . Assuming 
that the small grains are characterized by the same rate coef- 
ficient for momentum transfer through elastic scattering with 
neutrals as metal ions ((<rv) ± s=s 1.7xl0 _9 cm 3 s _1 ; see War- 
die & Ng 1999), we obtain rj w 0.07Cn%~ 1/4 (a/ 5 A )~ 1/4 in 
this density regime. Ambipolar diffusion will dominate Ohmic 
diffusivity so long as the grain Hall parameter (3± (the ratio 
of the cyclotron frequency to the collision frequency with neu- 
trals) exceeds 1 . Adopting a fiducial magnetic field strength at 
r = 10 AU of B = 0.3 G, we find 0± w 0.7(B/0.3G) /n i2 , which 
we expect to climb securely above 1 as the density continues to 
decrease with increasing r. We thus infer that both the low- and 
the high-density regimes of our model collapsing cores corre- 
spond to ambipolar diffusion that is amenable to self-similar 
scaling. The intermediate density regime will depart from this 
scaling (e.g., Nishi et al. 1991), and it is also likely to involve 
significant contributions from the Hall current term in Ohm's 
law (e.g., Wardle & Konigl 1993; Wardle & Ng 1999). Nev- 
ertheless, given that our bracketing values of the parameter r\ 
vary by only one order of magnitude (assuming that £n lies in 

the one used in CCK by a factor 
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the range <~ 1-10) across ~ 8 orders of magnitude in density, 
we may expect our diffusivity parameterization to yield quali- 
tatively meaningful results on all scales between <~ 10 AU and 
- 10 4 AU. 

Just as in CCK, we take the initial (f = 0) state of our model 
core to coincide with the endpoint of the pre-PMF collapse 
phase, which is characterized by a radial scaling of the sur- 
face density E and the vertical magnetic field component B z 
oc r~ x and by spatially constant radial inflow velocity V r (t = 
0) = ii()C (where C is the isothermal speed of sound) and mass 
accretion rate Mq = —A(uqC 3 /G). In the presence of rota- 
tion, this state also exhibits a uniform azimuthal speed. Since 
the latter is typically much smaller than the speed of sound 
C«0.19(r/10K) 1 / 2 kms- 1 , we normalize our initial state as 
in CCK based on the nonrotational collapse simulations of CK 
(mo ~ — 1 and A 3). We determine the initial value of V^/C 
(which we denote by vo) from the expression 



v 



Afl h C 

Vgb, 



0.15 



ref 



2xl0- 14 rads- 1 

/ Bref 



(2) 

^O.^kms- 1 J \30^G ' 

(see Basu 1997), where CI/, and B re f are, respectively, the uni- 
form background angular velocity and magnetic field strength. 
Based on the range of measured core angular velocities (see 
§ 1), vo could be a factor ~ 5 larger or smaller than the fiducial 
value adopted in equation (2). 

We assume that magnetic braking, which regulates the core 
rotation prior to its collapse, remains the dominant angular mo- 
mentum transport mechanism also during the subsequent evo- 
lution. [A centrifugally driven disk wind may, however, dom- 
inate the angular momentum transport in the rotationally sup- 
ported disk. In § 3 we argue that this is, indeed, a likely pos- 
sibility, and in Appendix C we show that this mechanism can 
be incorporated into our model without significantly modifying 
the basic formulation.] We adopt the approach articulated by 
Basu & Mouschovias (1994) for the pre-PMF collapse phase. 
The torque per unit mass on the slab-like core (of surface mass 
density S) can be approximated by rBJifo/lirE,, where B z is 
the vertical field component at the midplane and B^ jS is the az- 
imuthal field component at the surface z = H (assuming H <C r). 
The latter can be estimated from the relation 

where ty(r) is the magnetic flux threading the core within a 
radius r, V, \<p = ril, is the azimuthal speed of the core ions, 
and VA,ext is the speed with which the torsional Alfven waves 
that effect the magnetic braking propagate in the external 
medium. Equation (3) is the same as equation (26) in Basu & 
Mouschovias (1994), where VA,ext is also treated as a constant, 
except that we substituted the ion speed for the bulk (neutral) 
speed to take account of the possibility that the ions (into which 
the magnetic field is frozen) are not well coupled to the neu- 
tral core component (see Mouschovias & Paleologou 1986 and 
Konigl 1987). 

The assumption of a constant VX.ext is consistent with the- 
oretical models of magnetically supported clouds as well as 
with empirical data, which indicate that VA, ex t ~ lkms~' in 
molecular clouds over at least 4 orders of magnitude in den- 
sity (<~ 10 3 -10 7 cm~ 3 ; e.g., Crutcher 1999). This assumption 

3 As we discuss in § 3.2.3, we also do not expect the radial surface field component 



is required for the construction of a self-similar model. In 
order for the steady-state expression (3) to be applicable, the 
Alfven travel time in the external medium over the initial ra- 
dius of the cloud should be less than the evolutionary time 
scale <~ r/\V r \ (which, for a rotationally supported disk, is es- 
sentially the local magnetic braking time). Our representative 
rotationally supported disk models typically have |V r | < 0.1C 
(with \V r \ — > as r — > 0), whereas, for a nominal temperature of 
10K, V A ,exi/C w 5. These estimates indicate that the assump- 
tion of rapid braking should not lead to gross errors even if it 
is not everywhere strictly correct. In our models we use the pa- 
rameterization VA,ext = C/a, with a = 0.1 adopted as a typical 
value. 

We will be interested in evaluating the effects of magnetic 
braking in regions where rotation already plays a dynamically 
significant role, so that f2, » Qj,. We can therefore drop fij 
in equation (3), which will make it possible to incorporate this 
equation into a self-similar formulation. The anticipated in- 
crease in f2 ( - with decreasing r (oc r~ 3 /2 when a Keplerian disk 
is formed) will in general make the dominant field com- 
ponent at the surface as the central point mass is approached. 
(This situation is not encountered in the pre-PMF phase; see 
Basu & Mouschovias 1994.) We expect, however, that vari- 
ous magnetohydrodynamic instabilities (in particular, internal 
kinks) will intervene to prevent the azimuthal field component 
from greatly exceeding the poloidal components. 3 We therefore 
impose a cap on the azimuthal field in the form \B$. S \ ^ SB Z . In 
our models we usually set S = 1 — this choice, in fact, also 
corresponds to the typical value of \B^\/B Z obtained in our 
model for a rotationally supported, diffusive disk when (as ex- 
pected) the vertical angular momentum transport is dominated 
by a centrifugally driven wind (see § 4 and Appendix C). 

2.2. Basic Equations 

We represent the collapsing core as a thin disk surrounding a 
point mass M c . We use cylindrical coordinates (r, <f),z) centered 
on the point mass, with the disk midplane given by z = 0. The 
disk has a column density S, mass density p, and half-thickness 
H defined by E = 2Hp. The total mass M(r) enclosed within 
the radius r is then given by M(r) =M c + 2n E(r')r' dr' , and 
we denote its time derivative by M. The disk velocity field 
V(r) has radial and azimuthal components, with the latter giv- 
ing rise to a specific angular momentum J = rV^. We assume 
that the disk is threaded by an open magnetic field that is sym- 
metric with respect to the midplane [i.e., B r = at z = and 
B r , s {H) = —B r j(—H) at the disk surfaces (z = ±H), and simi- 
larly for B$\. The magnetic flux enclosed within the radius r is 
given by *(r) = * c + 2tt j^B z (r')r' dr' , where * e denotes the 
flux trapped inside the central point mass. 

In Appendix A we write down the basic disk equations, in- 
tegrate them over z while retaining all terms of order H/r, and 
then go on to justify the omission of some of the £)(H/r) terms 
in the interest of greater simplification. The result of these ma- 
nipulations is the following set of equations, 



«9£ 1 d 

ir + - 7r ('"V r S) = 
at r or 



(4) 



dV r „dV r C 2 <9£ B 7 

-+V r " =Sr 1 — 

dt dr Sr E 8r 2ttE 



B r , —H- 



dr 



J 2 

+ -, (5) 



dJ 
dt 



dJ rB 7 B 



' dr 



z D 4>,s 



2ttE 



(6) 



to greatly exceed the vertical field B z under realistic circumstances. 
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and 



EC 2 
~2H 



^2 , GM cP H 2 
2 GE + ^r^~ 



1 

8^ 



Or 



(7) 



which express mass, radial momentum, and angular momentum 
conservation, and vertical hydrostatic equilibrium, respectively. 

Under the assumption that the magnetic field evolution is 
governed by ambipolar diffusion, one can regard the magnetic 
field lines as being frozen into the ions. The latter move with 
velocity V, and drift with respect to the bulk (neutral) gas com- 
ponent at the drift velocity Vp = V, — V. By balancing the 
Lorentz force on the ions with the ion-neutral collisional drag 
force, one can solve for the components of V D : 



2ttS 
Tni B z B 



H 



dr 



2tt£ 



(8) 



(9) 



We then have 
~dt 

As in CCK, we adopt the monopole approximations for g r 
and B r „ 



-2nrV iir B z = -2nr (V r + Vb,r) B z 



(10) 



gr 



B rs 



GM(r,t) 
V 2 

*(r,0 



2irr 2 



(11) 



(12) 



which considerably simplify the calculations and are not ex- 
pected to introduce any significant errors (see also Li & Shu 
1997, Saigo & Hanawa 1998, CK, and Tsuribe 1999). 

For the azimuthal magnetic field component, we use equa- 
tion (3) in the limit V, ^ /r ^> fij together with equation (9), and 
impose the bounding condition on \B^\ discussed at the end of 
§ 2.1. We then obtain 



Ba 



mm 



1 + 



*B 7 



7rr 2 V A ,ext V 27r 2 r 2 £V A ,ext 



:5B 7 



(13) 



2.3. Self-Similar Equations in Nondimensional Form 

We introduce a similarity variable x and a set of nondimen- 
sional flow quantities that depend only on x: 

x = r/Ct, (14) 

H(r,t)=Cth(x), E(r,f) = (C/2nGt ) a(x) , (15) 

V r {r,t) = Cu{x) , V^{r,t) = Cv(x) , (16) 



g T {r,t) = (C/t)g(x) , J(r,t) = C 2 tj(x) , 



(17) 



M(r,t) = (C 3 t/G)m(x), M(r,t) = {C 3 /G)m{x) , (18) 
B(r,t) = (C/G 1/2 f)b(x),*(r,f) = {2irC 3 t /G l/2 )ip(x) . (19) 

These definitions extend the set used in CCK to the rotational 
case (which involves the dimensionless azimuthal speed v and 
specific angular momentum j) — note, however, that the vari- 
ables h, u, a, and /i = a Jb z in the present paper were denoted 
by h, v, a, and A, respectively, in CCK. To help relate our re- 
sults to real protostellar systems, we observe that, for our fidu- 
cial value of C (= 0.19kms _1 ), x = 1 corresponds to a distance 



r :=y 6 x 10 15 cm (= 400 AU) when t = 10 4 yr (the characteris- 
tic age of a Class YSO) and to a distance r«6x 10 16 cm 
(= 4000 AU) when t = 10 5 yr (the characteristic age of a Class 
I YSO). 

CCK demonstrated that ambipolar diffusion can be incorpo- 
rated into a self-similar model when p, oc p 1 / 2 . When rotation 
is present, we find that magnetic braking can similarly be in- 
cluded if VA,ext = const. The nondimensional model parameters 
that control the strength of these two effects are r\ and a, re- 
spectively (see § 2.1). [A self-similar formulation is also pos- 
sible if angular momentum transport is due to an "a viscosity," 
v = a ss CH (Shakura & Sunyaev 1973), when C and the param- 
eter ass (rS 1) are constant. This was noted by Tsuribe (1999) 
in the special case where H has the value appropriate to a Ke- 
plerian disk, and by Mineshige & Umemura (1997) in the case 
when H has the value appropriate to a self-gravitating disk (al- 
though the latter authors required assH/r, rather than ass, to 
be a constant).] With the above expressions for the similar- 
ity variables, and defining also w = x — u for convenience, the 
structure equations can be rewritten in the following nondimen- 
sional, self-similar form: 



du 
dx 



1 da 1 

a dx x 



, 1 da b z . , 

(l-W 2 )-—=g + -[ b r . 

v ' a dx a 



db z 
dx 



■2 2 
J w 



+ ^T + 



(20) 



(21) 



ip - xwb z + rjxb 2 h l/2 a- y2 [b^ - h^j = , (22) 



dj_ 

dx 



w 1 (j-xb z b<f,, s /a) , 



(23) 



. [2aV7'/ 1 2ar l h 1 / 2 ipb z \ ' „ l 
b^ = -mm I — p (l + J z - ) ; Sb z } , (24) 



a 3/2 x 2 



m = xwa , 
m = —xua , 



dm 
dx 

dx 



■ xa 



= xb z , 



-m/x 2 



and 



b r . s = ip /x 2 . 



(25) 
(26) 

(27) 

(28) 
(29) 

(30) 



The disk half-thickness h can be found from the vertical hy- 
drostatic equilibrium condition, which yields 

(^T " h > '^) # + K + b l, + « 2 )h-2a = 0. (31) 
The solution of this quadratic equation is 



h = 



1+1 + 



6m c 



1/2" 



(32) 



2m c 

where m c = m c — x i b r . s (db z /dx)/a and a = a+ (b 2 s + b^ ^/a. 
Equation (32) implies that h — > 2/a, (2/m e ) 1//2 x 3 / 2 , and 



2a/[b 2 s + b 2 bs ] in the limits where self-gravity, central-mass 
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gravity, and magnetic stresses, respectively, dominate the verti- 
cal squeezing of the disk. 4 

As discussed in § 2.1, the initial it = 0) conditions [which 
in the self-similar model also represent the outer asymptotic 
(r — > oo) values] correspond to a collapsing core just before 
PMF. Thus we require 

A a 

a ^ — , — » — , u — > uo, v — > vo as x — > oo , (33) 
x ^ 

where A, /Uo, Wo, an d v o are constants. As discussed in § 2.1, 
our fiducial values for A (= 3) and mo (= —1) are the same as 
those used by CCK; they are based on the results of numerical 
simulations of nonrotating cores and are compatible with ob- 
servations. We similarly adopt the CCK choice fio = 2.9. Our 
reference values for the uniform initial rotation vo are based on 
the range of measured core angular velocities and typical cloud 
parameters that enter into the estimate (2). 



Our "standard" solutions, corresponding to the fiducial pa- 
rameter values, are presented in § 3.2.1 for the IMHD case and 
in § 3.3.1 for the AD case. We also present solutions in the 
limit when the initial rotation is very fast and magnetic braking 
is either nonexistent (§ 3. 1) or comparatively inefficient (§ 3.2.2 
and § 3.3.2 for the IMHD and AD cases, respectively). Under 
these conditions the rotationally supported disk becomes fairly 
extended and its mass greatly exceeds that of the central object. 
Another "extreme" type of a solution is obtained in the oppo- 
site limit, when the initial rotation is braked so efficiently that 
j is effectively reduced to zero at a finite value (Xj) of x. In this 
limit, too, we present solutions for both the IMHD (§ 3.2.3) and 
AD (§ 3.3.3) cases. These configurations are effectively nonro- 
tating (and thus resemble the CCK solutions) for x < xj, and 
they do not feature either a centrifugal radius or a Keplerian ac- 
cretion disk; instead, matter accretes directly onto the central 
object at a supersonic speed. 



3. RESULTS 

We have studied three distinct cases of rotational core col- 
lapse, which we present in this section in order of increasing 
complexity. They are: (1) A purely hydrodynamical collapse, 
with no mechanism for angular momentum transport (§ 3.1); 
(2) An ideal-MHD (IMHD) collapse, which incorporates mag- 
netic braking but does not include any magnetic field diffusivity 
that could prevent the buildup of a central magnetic monopole 
(§ 3.2); (3) an MHD collapse that includes ambipolar diffusion 
(AD; § 3.3). The latter model is able to reproduce the basic ob- 
servational features of rotationally supported protostellar disks 
and their central YSOs. 

Our models have many features in common, arising from the 
basic interplay between the centrifugal force and gravity. Ini- 
tially and at large distances from the center (i.e., for large values 
of x), gravity is stronger than the centrifugal force and matter 
falls in at a high, supersonic speed. In this regime m ^> m c , 
and self-gravity dominates over attraction by the central ob- 
ject. A very different behavior characterizes small values of x. 
Typically, the innermost region constitutes a Keplerian accre- 
tion disk, where gravity approximately balances the centrifu- 
gal force and where the infall speed is low and subsonic. The 
transition between the supersonic and subsonic inflow regimes 
is achieved through a shock, which is located at the point (x c ) 
where the infalling matter encounters a centrifugal barrier. In- 
ward of this shock the column density increases significantly, 
forming a dense accretion disk. The location of the centrifu- 
gal shock roughly coincides with the centrifugal radius, which 
is the largest radius where the gravitational (— m/x 2 ) and cen- 
trifugal ij 2 /x 3 ) accelerations are in approximate balance (e.g., 
Basu 1998). In fact, we find that the point where the equal- 
ity j 2 (x) jm(x) — x is first satisfied as the matter flows in yields 
a value that is only a slight (~ 4%) overestimate of the actual 
shock location in both the IMHD and AD cases when (as is the 
case in our typical models) m(x c ) w m c . The value of the cen- 
trifugal radius is sometimes estimated from the t = expres- 
sions for the surface density and angular momentum profiles. 
This estimate could be rather inaccurate since it does not ac- 
count for the mass accumulation at the center or for the loss of 
angular momentum that the infalling matter experiences even 
before it reaches the centrifugal radius. We discuss more accu- 
rate estimates of x c in the following subsections. 

4 CCK employed an expression similar to equation (31), except that they omitted 
but with m c and a replaced by m c and a, respectively. 



3.1. Nonmagnetic Rotational Collapse 

This problem was already treated in detail by Saigo & 
Hanawa (1998); it is included here for reference and compari- 
son, and to serve as a starting point. 

In the absence of magnetic fields or of any other means 
of angular momentum transport, the self-similar collapse is 
characterized by j/m = <&, where $ is a constant fixed by 
the initial conditions (<£> = vq/A). The inner asymptotic so- 
lution (x — ► 0) is given by a = E/x and m = Ex, with E 

I + (1 - 4$ 2 ) 1 / 2 ]/(2$ 2 ). Numerically, we proceed by ap- 
plying the initial conditions at x max = 10 2 and integrating the 
differential equations toward smaller values of x until the ex- 
pected location of x c (see eq. [36] below) is reached. There we 
impose the isothermal shock jump condition described in Ap- 
pendix B.l, and continue the integration until the asymptotic 
solution is approximately fulfilled. The exact value of x c is ob- 
tained iteratively by enforcing the matching to this asymptotic 
solution to a high precision. 

Figures 1 and 2 show profiles of mass, column density, and 
radial speed for two kinds of initial rotation: a typical (vo = 0.1) 
and a rapid (vo = 1) one. In both cases the centrifugal shock is 
very abrupt, producing a large increase in column density and a 
sudden decrease in inflow speed. In the case with faster initial 
rotation, the effect of the centrifugal barrier is so strong that it 
forces the radial velocity u to change sign below x c to u > 
(indicating a backflow), as shown in Figure 2. The centrifu- 
gal shock thus separates regions of infall and backflow; it stays 
stationary in self-similar coordinates, implying that it moves 
outward in physical space. This phenomenon was already dis- 
cussed in Saigo & Hanawa (1998). By trying a few other values 
of vo (but keeping the other parameter values unchanged), we 
have found that vo = 0.5 is too small to produce a backflow but 
that vo = 0.8 is already large enough. 

An interesting behavior exhibited by these solutions (and 
standing out particularly in Fig. 1) is the development of a mass 
"plateau" when the infall speed starts to increase above its outer 
asymptotic (x ^> 1) value. This behavior is generic to all ex- 
cept our fast-rotation solutions (§ 3.2.2 and § 3.3.2) and can be 
understood as follows. According to equations (25) and (26), 
m = ax(x — u) and m = —oxu. These equations are consistent 
with the outer asymptotic solution m=Ax = ax 2 and u = uq < 
(eq. [33]) only so long asi> \u\ (in which case m 3> m). This 

magnetic terms. Their solution for h thus has the same form as equation (32), 
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FIG. 1. — Self-similar solution for a nonmagnetic rotational collapse. The 
variations of the normalized radial infall speed —u, surface density a, and total 
mass m are plotted as functions of the similarity variable x. The initial (or outer 
asymptotic boundary) conditions correspond to the parameter values vo = 0.1, 
A = 3, and uq = —1. The centrifugal shock is located at x c = 7.5 X 1()~ 3 . No 
central mass forms in this case on account of the assumed lack of any angular- 
momentum transport mechanism. 
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FIG. 2. — Same as Fig. 1, except that the azimuthal-velocity parameter v is 
increased from 0.1 to 1. Note the change in the sign of u across the centrifugal 
shock in this case and the comparatively large value of x c (=0.63). 



m p i W (A/x p i)Xpi(Xpi - M ) 

> x, and hence m 



inequality breaks down when x becomes small enough that it 
drops below \u\. This point marks the approximate outer edge 
of the mass plateau, and we label it x p i. Its value can be esti- 
mated by setting 

*pi~|«o|, (34) 

although typically \u\ is already somewhat larger than \uq\ at 
that point. The plateau mass can then be approximated by 

2\uo\A. (35) 

For x <C x p i, \u\ » x, and hence m w m and (using eq. [27]) 
dm/dx = (a < m/x The latter inequality accounts for the ap- 
pearance of the plateau in the logarithmic plots of m vs. x. If 
x c /x p i is sufficiently small, then, in the plateau regime, \u\ ap- 
proaches the free-fall value corresponding tomw m p \. Hence, 
once \u\ comes to exceed x, their ratio increases rapidly with 
decreasing x, and the mass plateau is quickly established. The 
derivation of equations (34) and (35) is independent of the de- 
gree of magnetization of the flow, and the expression (35) there- 
fore also provides an excellent first guess of m c in our fiducial 
magnetic solutions that involve point-mass formation. These 
results do not, however, apply to the fast-rotation cases con- 
sidered in § 3.2.2 and § 3.3.2, where a plateau does not form 
because the centrifugal shock is established before x becomes 
<C \u\, and where the central mass is typically smaller than that 
of the disk. 

Setting j/m = $ and m s=a m v \ in the expression x c 
(j/m) 2 m, we infer 

x c tt2\u \vl/A. (36) 

We have used this expression to obtain the initial guess in our 
numerical calculation. For the parameters of Figure 1 , this gives 
x e m 7 x 10~ 3 , which is seen to be an excellent approximation 
to the actual value. This estimate is only applicable if x c < x p i, 
i.e., for vo ^ 0V2) '^ 2 - This condition is still marginally satis- 
fied in the solution presented in Figure 2. 

The inner asymptotic solution of this model precludes the 
formation of a point mass. As the azimuthal speed v cannot di- 
verge as x — > 0, j = vx must vanish at the origin. To form a cen- 
tral object with m c ^ 0, a braking mechanism that allows j/m 
to decrease as x — > is required. The absence of such a mech- 
anism explains the linear decline of m with x at small values 
of x (which for the parameters used in Figs. 1 and 2 effectively 
starts already at x c ). It may also be expected that a braking 
mechanism will lead to a more gradual radial deceleration, re- 
ducing or even altogether preventing the backflow exhibited by 
the rapid-rotation solution shown in Figure 2. In the following 
subsections we examine to what extent these expectations are 
met in the presence of magnetic braking. 

3.2. Ideal-MHD Rotational Collapse 

In this section we present an IMHD model for magnetic brak- 
ing based on the formalism described by equations (20)-(31), 
in which we set 77 = 0. 

The existence of a braking mechanism leads to a qualitative 
change in the character of the solutions. In contrast to the non- 
magnetic solution presented in § 3.1, j/m is not a constant in 
this case, and magnetic braking allows a central mass to form 
despite the initial presence of rotation. 

Flux freezing under the ideal-MHD assumption implies that 
b z = a/fi and ip — m//i, where \i, the nondimensional mass-to- 
flux ratio, is a constant (equal to its initial value fio). It is also 
convenientto introduce the auxiliary constant p,= (1 — /i -2 ) 1 / 4 , 
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which is useful for describing the effect of the magnetic "dilu- 
tion" of the gravitational field (Shu & Li 1997). 

There are two qualitatively different sets of asymptotic re- 
lations for the inner region (x — > 0) of this model. The first 
set, presented in equations (37)-(42) below, applies to the con- 
figurations described in § 3.2.1 and § 3.2.2. The second set, 
presented in § 3.2.3, applies only in the strong-braking case. 
Both of these sets describe the formation of a dynamically dom- 
inant central mass and a split-monopole field. However, the first 
set represents a rotationally supported disk, whereas the second 
one depicts a nonrotating, supersonic inflow. 

The disk-like inner asymptotic solution is 



m = m = 


m c , 


(37) 


j = 


ftWJ 2 ^ 2 , 


(38) 


— U = w = 


(2«5) 1 / 2 myV- 1 M-^ 1/4 , 


(39) 


= b z = u/\x = 




(40) 


b r . s = iP/x 2 = 


m c /(iix 2 ), 


(41) 


h = 2n/b% = 


M M V 5 /4 (2/(5) i/2 x ii/4. 


(42) 



These equations represent an essentially Keplerian disk; in par- 
ticular, the expression for j differs from that of a purely Kep- 
lerian system only by the magnetic-dilution factor p 2 [which, 
however, remains close to 1 (« 0.94) for our fiducial value of 
|Uo]. The j profile corresponding to such a disk implies that 
j/m decreases with decreasing x. The magnetic braking re- 
quired to account for this change from the initial distribution 
fixes the strength and radial power-law index of the magnetic 
field components b z and b^, s . The flux-freezing condition then 
determines the value of a, and the mass conservation relation, 
in turn, yields the expression for u. 

The behavior of the IMHD solutions is determined primar- 
ily by the magnitude of the initial rotation and the strength of 
the magnetic braking. We consider three qualitatively different 
cases: a typical collapse in which both effects are relevant over 
most scales (§ 3.2.1); a limiting case corresponding to weak 
braking and fast rotation (§ 3.2.2); and another limiting case in 
which very strong braking is capable of eliminating rotation al- 
together by the time the flow reaches a finite distance from the 
origin, resulting in the formation of an inner nonrotating-infiow 
zone (§ 3.2.3). 

The method used to solve the differential equations is essen- 
tially the same as that used in § 3.1. For the initial guess of 
the value of x c we again use the estimate (36). Although mag- 
netic braking does operate in this case, the deviation of j/m 
from its x ^> 1 asymptotic value remains relatively small in the 
range x > x c , and the approximation obtained for the nonmag- 
netic case remains fairly good (as evidenced by the fact that the 
value of x c in the solution presented in Fig. 3 differs by a fac- 
tor of less than 2 from its value in Fig. 1). The discontinuities 
in the variables a, b z , and u across the shock are estimated us- 
ing the appropriate shock jump conditions. For the typical case 
(§ 3.2.1) the shock occurs in a region where h m 2a/b 2 s and the 
"magnetically squeezed shock" jump conditions described in 
Appendix B.3 are applicable. In the fast-rotation case (§ 3.2.2) 
the shock is located sufficiently far away from the origin that 
a generalized (to include magnetic effects) isothermal shock 
(Appendix B.2) provides a better approximation. (The issue of 
shock jump conditions is not relevant to the strong-braking in- 
flow considered in § 3.2.3, since no centrifugal shock forms in 
that case.) As in § 3.1, a successful convergence of the solution 
at x < x c to its asymptotic (x — > 0) form requires a correct se- 



lection of the point x c , which is achieved iteratively. (Incorrect 
choices of x c produce numerical solutions with a spontaneous 
singularity at some point < x < x c , similar to those found in 
Li 1998.) As the value of the central mass is not known before 
the numerical integrations start, we adopt m c w m p i (eq. [35]) 
as an initial guess that is refined iteratively until convergence is 
attained. 

3 .2. 1 . Fiducial Solution 

The results for this case are presented in Figure 3. For the 
fiducial parameters vo = 0.1, a = 0.1, and S = 1, rotation and 
braking are both moderate. Magnetic braking is strong enough 
to prevent shock formation in the region where self-gravity 
dominates (as in § 3.2.2), but it is not so strong as to prevent 
the shock from forming altogether (as in § 3.2.3). The centrifu- 
gal shock thus occurs well inside the region where the central 
mass dominates the gravitational field, with m(x c ) « m c . The 
disk that forms for x < x c is thus essentially Keplerian and has 
an angular momentum profile j oc x x l 2 . The large inclination of 
the magnetic field lines to the disk surface indicated by the in- 
ner asymptotic solution (b ns /b z oc x~ 3 / 4 as x — > 0; see eqs. [40] 
and [41]) suggests that the disk would likely drive a centrifu- 
gal wind from its surfaces, which would reduce m c . We neglect 
this effect in the present discussion, but we consider it in § 4 
and Appendix C in relation to the ambipolar-diffusion fiducial 
solution. 

The evolution of the inflow can be characterized by dividing 
the range of x into four regions, which are, in order of decreas- 
ing distance from the center: 

Self-gravity-dominated region. For the largest values of x, 
the initial value of the ratio j/m (— 0.033) is still conserved, 
and m, which is still much larger than its central value, is pro- 
portional to x. 

Free-fall region. After a transition region, where j/m de- 
creases slightly from its initial value on account of magnetic 
braking, both m and j enter a plateau (at x w x v \\ see eq. [34]). 
They become independent of x, with values m v \ = 6.2 and 
j'pi =0.18 (corresponding to j/m = 0.028). These computed 
values can be substituted into the expression x c m j P i 2 / m pi to 
provide a refinement of the estimate (36) for the location of 
the centrifugal shock. As can be seen from a comparison of 
this value with the computed result (listed in the caption to Fig. 
3), the approximation is excellent (as is the correspondence be- 
tween the value of m p i and that of m c ). Near the inner edge of 
this region, the centrifugal force starts to exceed gravity. 

The centrifugal shock. The outer edge of this region is 
where the variables a, u, and b z undergo the discontinuities 
prescribed by the shock jump conditions of Appendix B.3. 
Just inside the discontinuity the column density becomes very 
large, but it then rapidly decreases toward its asymptotic near- 
Keplerian value. The infall speed decreases rapidly, but not so 
strongly as to produce any backflow (i.e., u remains < for 
all x). It is in this postshock transition zone that the angular 
momentum makes the final adjustment to its asymptotic, mag- 
netically diluted, Keplerian value: this sudden bout of braking 
is triggered by the spike in b z that results from the shock com- 
pression. The column density spike that occurs in this region 
is associated with a small (but finite) mass, Am w 0.02, which 
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FIG. 3. — Fiducial solution for an ideal-MHD rotational collapse. The vari- 
ations of the normalized radial infall speed — u, surface density a, total mass 
m, mass accretion rate m, specific angular momentum j, and z component of 
the magnetic field b, are plotted as functions of the similarity variable x. The 
initial (or outer asymptotic boundary) conditions correspond to the parameter 
values v = 0.1, A = 3, u = — 1, and ^ = 2.9, whereas the magnetic-braking 
model parameters are a = 0.1 and 5=1. The centrifugal shock is located at 
x c = 4.9x 10~ 3 , and the horizontal scale in its vicinity has been expanded to 
show details of the postshock transition zone. The central mass is m c = 6.0. 



x = r/Ct 

FIG. 4. — Fast-rotation solution for an ideal-MHD rotational collapse. The 
variations of the normalized radial infall speed — u, surface density a, total 
mass m, and specific angular momentum j are plotted as functions of the sim- 
ilarity variable x. The model parameter values are the same as in the fidu- 
cial case (Fig. 3), except that the azimuthal-velocity parameter i> is increased 
from 0.1 to 1.5. In this case the centrifugal-shock radius is comparatively large 
(x c = 1.48) and the central mass is rather small (m c = 0.57). A backflow layer 
(u > 0) is present just inside the centrifugal shock: within this layer, the veloc- 
ity curve depicts +u instead of — u. 



could be pictured as a massive ring located at the edge of the 
Keplerian disk. After undergoing fast variations as it transits 
through this region, the flow merges seamlessly into the next, 
smoothly varying zone. 

Magnetically diluted Keplerian disk. In this region, the 
asymptotic equations (37)-(42) are satisfied. The rotationally 
supported disk has a mass nid s=a 3.2% of m c and is in an al- 
most perfect dynamical equilibrium, as the inflow speed is low 
(<0.2C). 

3.2.2. Fast Rotation 

If the initial azimuthal speed vn is high and the braking pa- 
rameters a and 8 are moderate or small, then the value of j can 
become quite high, which can lead to a large value of the cen- 
trifugal radius x c . In an extreme case, the centrifugal shock may 
occur so far away from the center that it is located within the 
self-gravity-dominated region, where m^>m c . 

Figures 4 and 5 show two solutions, corresponding to a fast 
initial rotation (vn = 1.5) and two different values of the braking 
parameter a (0.01 and 0.1, respectively; 6 = 1 in both cases). 
The centrifugal shock in these solutions is so strong that the ra- 
dial speed u changes sign across the shock. This is similar to 
the backflow observed for large values of j in the absence of 
braking (see Fig. 2) and is consistent with the expectation that 
the braking must be strong enough to qualitatively change the 
behavior of a rapidly rotating inflow. The backflow region in 
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X = r/Ct 

FIG. 5. — Same as Fig. 4, except that the braking parameter a is de- 
creased from 0.1 to 0.01. This results in a wider backflow layer and a smaller 
Keplerian-rotation region. The centrifugal shock location is not much changed 
{x c = 1.53), but the central mass is strongly reduced (m c = 0.05). 
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this case is clearly seen to occupy only a finite range in x and 
not to reach all the way to the origin: it represents an annulus, 
interior to which the infall resumes. 

The large values of x c result in extended, and hence fairly 
massive, disks. On the other hand, the comparatively weak 
braking inhibits mass accumulation at the center: the central 
masses obtained in these solutions (listed in the captions to 
Figs. 4 and 5) are much smaller than the mass derived for the 
fiducial parameters in § 3.2.1 (which, in turn, is comparable to 
m c in the nonrotating collapse model of CCK). The outer re- 
gions of the disks are therefore manifestly non-Keplerian. (The 
specific angular momentum approaches to within ~ 20% of 
the asymptotic, diluted-Keplerian value given by eq. [38] only 
when x decreases below <~ 0.1 \x c for the solution in Fig. 4 and 
below <~ 0.04x c for the solution in Fig. 5.) As can be seen from 
a comparison of Figures 4 and 5, a larger value of the braking 
parameter a leads to a reduction in j and results in the value 
of x c becoming smaller and that of m c larger. Correspondingly, 
the width of the backflow region is reduced and the Keplerian 
regime extends further out. These trends continue as the effect 
of rotation is diminished: by the time the parameters attain their 
fiducial values, the backflow region has completely disappeared 
and the inner asymptotic solution starts very close to x c . 

In summary, the collapse of fast-rotating cores that lack an 
exceptionally strong braking mechanism results in fairly mas- 
sive disks and in comparatively low-mass central objects. The 
disks are rotationally supported and largely self-gravitating, 
with only their innermost regions exhibiting Keplerian behav- 
ior. 

3.2.3. Strong Braking 

If braking is very efficient then the innermost region of the 
inflow can have effectively zero angular momentum even if the 
collapsing core starts out rotating. This situation is illustrated 
in Figure 6, which depicts a solution with a fairly respectable 
initial rotation (vo = 1) but with extreme values of the braking 
parameters (a = 5 = 10), which give rise to a strong surface az- 
imuthal field (see eq. [24]). In this case the angular momentum 
is reduced effectively to zero at xj w 0.2, before a centrifugal 
shock could occur, and thus no rotationally supported disk is 
formed. For x < xj, j = b^^ = 0, and as x — > 0, the inflow tends 
to the asymptotic form of the nonrotating magnetic solution de- 
rived by CCK (see their §3.3). In this limit the flow becomes a 
supersonic, magnetically diluted free fall onto the central mass: 

m = m = m c , (43) 

-u=[2m c {l-n- 2 )/x\ L ' 2 , 

<r = /j,b z = [m c /2(l — n~ 2 )x\ l l 2 , 

b rtS = ip/x 2 = m c j (fix 2 ) , 

(.2 _ ,.2r 



h = 2a/b 2 rs = y\2lm\(\ - ^" 2 )] 1/2 x' 'I 2 . 



(44) 
(45) 
(46) 
(47) 



The numerical solution converges to these asymptotic values 
already quite close to Xj (atx <~ 0.1). 

Equation (47) for the disk scale height shows the strong effect 
of magnetic squeezing by the radial magnetic field component 
in the asymptotic solution. (Since the radial surface magnetic 
field scales as x~ 2 , it dominates the magnetic terms in both the 
vertical and the radial force-balance equations.) The depen- 
dence of h on x (oc x 7 / 2 ) is much stronger than in the solution 
presented by CCK (h oc x 3 /2 ), in which the effect of magnetic 
pinching was neglected. However, apart from this difference, 
these two solutions are very similar; in particular, the value of 
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FIG. 6. — Strong-braking solution for an ideal-MHD rotational collapse. 
The variations of the normalized radial infall speed — w, surface density a, to- 
tal mass m, and specific angular momentum j are plotted as functions of the 
similarity variable x. Moderately fast initial rotation (v = 1) and large brak- 
ing parameters (a = 5 = 10) are assumed; the other parameter values are the 
same as in the fiducial solution of Fig. 3. In this case j decreases to a very 
small value at a finite distance from the origin (xj ~ 0.2), and no centrifugal 
shock is established. The central mass is m c = 6.06. 



the central mass quoted by CCK (m c = 6.1) is very close to 
the value that we derive (m c = 6.06). We note, however, that 
the extremely rapid decrease of h with radius that is implied by 
this solution is unlikely to be sustained in reality. For example, 
when the disk becomes sufficiently thin, the radial field compo- 
nents on opposite sides of the disk might reconnect, resulting in 
an upper bound on b v . Furthermore, the decrease in h will be 
arrested once the column density grows to a value for which the 
disk becomes opaque and the isothermal approximation breaks 
down. We do not consider these caveats in more detail since 
(as we point out at the beginning of § 3.3) we expect radial 
field diffusion to intervene (possibly even before the above ef- 
fects come into play) and prevent the IMHD configuration from 
being set up near the center. 

3.3. Ambipolar Diffusion-Dominated Rotational Collapse 

Using the x — > asymptotic relations for the IMHD solu- 
tion given by equations (37)-(42), it can be readily verified that 
the nondiffusive terms in the induction equation (22) are oc x° 
whereas the diffusive term (associated with the magnetic ten- 
sion force) is oc x~ 1//4 , so that the latter will become dominant 
near the origin. It follows that the IMHD solution cannot apply 
in the vicinity of the center in cases in which an accretion disk 
is predicted to form. Ambipolar diffusion must therefore be in- 
corporated into the model, and in this subsection we describe 
how the nature of the solutions changes when this is done. 

In the case of the asymptotic solution for a strongly braked 
rotational infall or a nonrotational collapse (eqs. [43]-[47]), the 
nondiffusive terms in equation (22) are still oc x°, but the dif- 
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fusive term scales as x 1 / 2 and remains subdominant as x — > 0. 
The asymptotic IMHD solution is thus self-consistent in this 
case, as we verified by checking that, in the model presented 
in Figure 6, the diffusive term is everywhere smaller than 10% 
of the fiux-advection term (assuming r\ < 1). This conclusion 
differs from that of CCK, a discrepancy that can be traced to 
their neglect of the magnetic squeezing term in the equation for 
h. However, as we noted in § 3.2.3, real disks are unlikely to 
become as extremely thin as the IMHD solution that incorpo- 
rates this term would imply. Furthermore, the accumulation of 
magnetic flux at the center, which is indicated by the IMHD 
split-monopole solution, would lead to a severe magnetic flux 
problem (see § 1) if allowed to proceed and is likely prevented 
from occurring by the action of Ohmic diffusivity in the disk or 
the YSO (see Li & McKee 1996). 

The ascendancy of the ambipolar diffusion term at low values 
of x means that ions and neutrals behave differently in that re- 
gion. In contrast, the initial (x ^> 1) conditions involve no field 
diffusivity and correspond to negligible drift between these two 
particle species. The study of these two regimes can be facil- 
itated by observing that, in most cases, the strong inequality 
b r .s 2> h(db z /dx) is satisfied for almost all values of x: the gra- 
dient of b z is never very large (except in shocks), and the disk is 
usually very thin. Under these conditions, one can treat the flux 
conservation relation (eq. [22]) as a quadratic equation for b z , 

r]ipx- 1 h 1/2 a- 3/2 b 2 -x{x-u)b z + ip = , (48) 

where we regard the variables ip, h, a, and u as given. Equa- 
tion (48) has two real roots, b z \ ow and b z ,hi g h, and we have found 
that in most of the solutions that we constructed they are well 
separated (b Z: \ 0VI <C £> z .Mgh)- This property makes it possible to 
obtain simple expressions for b z in the IMHD (large x) and 
AD (small x) regimes. Specifically, the lower root, which is 
applicable in the IMHD regime, can be evaluated by disregard- 
ing the quadratic term in equation (48): 



b z \ovi ~ ip/xw — aip/m. 



(49) 



The vertical field component in the AD regime corresponds to 
the larger root and can be approximated by omitting the con- 
stant term in equation (48): 

bz.high ~ ?T 1 x(er//z) 1/2 (m/'(/>) . 



(50) 



The transition between the IMHD and AD regimes can take 
place either gradually or abruptly. The latter case occurs when 
the flow is super-fast-magnetosonic and the flow variables are 
changed by passing through a shock. The occurrence of such an 
ambipolar diffusion-mediated shock was originally suggested 
by Li & McKee (1996) and subsequently verified in the nonro- 
tating collapse calculations of CK and CCK, which succeeded 
in resolving the shock structure. We generalize these results to 
the rotating-collapse case in § 3.3.1. The most striking feature 
of this shock is the abrupt jump in b z between the two roots 
(49) and (50) (the other flow variables are also affected, but 
less strongly). The jump represents a large local gradient in b z , 
which implies that the term h(db z /dx) in equation (22) can no 
longer be disregarded. In fact, it is this term that mediates the 
transition between the two roots of equation (48) and allows 
the shock structure to be resolved. The transition, however, re- 
mains very sharp, and equation (48) remains applicable outside 
a narrow neighborhood of the shock position x a . 

As in the IMHD solution presented in § 3.2, we find it con- 
venient to distinguish between the fiducial case (§ 3.3.1) and 
the limiting cases of fast rotation (§ 3.3.2) and strong braking 



(§ 3.3.3). In the fiducial case the braking is strong enough to 
place the centrifugal shock well within the AD shock (x c <C x a ). 
The centrifugal shock again disappears in the strong-braking 
case, but an AD shock is still present, as in the nonrotating, dif- 
fusive solution obtained by CCK. However, in the fast-rotation 
case the centrifugal shock is established when the flow is still in 
the IMHD regime, and hence for x <x c the inflow speed is low. 
In this case the transition between the IMHD and AD regimes 
occurs smoothly instead of in a shock (with the two roots of eq. 
[48] remaining close to each other), although it is still mediated 
by the h{db z /dx) term in the induction equation. 

To solve the self-similar equations as a boundary-value prob- 
lem, it is necessary to know the asymptotic behavior for x — > 0. 
For a flow that continues to rotate all the way down to the ori- 
gin (as in the fiducial and fast-rotation solutions), the behavior 
in this limit is given by 

m = m = m c , (51) 

(52) 
(53) 



m\/ 2 x 1 ' 2 



■ w = {m c /ai)x l/2 , 
(j = j2 V /35)(2m c y/ 2 



-3/2 



(J\X 



b r . s 



-(277/35)" 2 ] 1 / 2 

-3/2 

5 

-W* = [«c 3/4 /(2*) 1/2 >- S/4 , 
= ^/x 2 = (4/3)b zl 

h = {2/[l + (2 V /38) 2 ]m c } l / 2 x^ 2 . 



(54) 
(55) 
(56) 
(57) 



This solution represents a Keplerian disk in which all the field 
components have comparable magnitudes. In contrast with 
the corresponding IMHD solution with its split-monopole field 
configuration (see eqs. [37]-[42]), in this case the field-matter 
decoupling brought about by the ambipolar diffusion results in 
a field that is not strong enough to either dilute the Keplerian 
rotation or to dominate the vertical compression. [The vertical 
squeezing of the disk in the asymptotic AD solution is dom- 
inated by the tidal and self-gravity forces, which contribute 
in the proportion 1 : (2rj/3S) 2 .] The asymptotic solution in 
the strong-braking limit, which (as in the corresponding IMHD 
case) involves a flow that is nonrotating near the origin, is given 
in §3.3.3. 

The solutions presented in this subsection have been obtained 
by using the following numerical procedure: 

1 . Solve a simplified system of equations by disregarding 
the db z /dx terms in equations (21), (22), and (31) to ob- 
tain approximate values for the five variables a, u, ip, j, 
and b z at a selected matching point x m . We have found 
that, to assure numerical stability, it is best to choose this 
matching point in the region where the effects of am- 
bipolar diffusion are just starting to become dominant. 

2. Integrate the full system of equations both outward and 
inward from the matching point. In the outward direc- 
tion, the solution reaches out to the point x max = 10 2 , 
where it must match the prescribed initial conditions of 
the flow (given by the parameters A, uq, m>, and Vo). 
In the inward direction, the flow must match the x — > 
asymptotic solution (given by eqs. [51]— [57] in the fidu- 
cial and fast-rotation cases, and by eqs. [66]-[71] in the 
strong-braking case), which also determines the value of 
m c . Except in the strong -braking limit, the flow must 
also incorporate a centrifugal shock (which is located 
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inside the matching point for the fiducial solution and 
outside it in the fast-rotation case). The discontinuities 
that some of the flow variables experience at x c are han- 
dled using the appropriate shock jump conditions (Ap- 
pendix B). 

3. Iterate on the values of the five flow variables at the 
matching point x m as well as on those of the central mass 
(m c ) and (if applicable) of the centrifugal-shock position 
(x c ) until convergence is reached. 

3.3.1. Fiducial Solution 

Our "standard" solution, corresponding to rj = 1, vo = 0.73, 
a = 0.08, and 5 = 1, is presented in Figure 7. In this case 
the initial rotation is not very fast and the braking is moderate 
(see § 3.2. 1), so the ambipolar-diffusion shock is located further 
away from the center than the centrifugal shock (x a /x c w 30). 
To more clearly separate the different flow regimes, we also 
show a solution where vo is reduced to 0.18 (Fig. 8a), which 
has the effect of reducing x c by a factor ~ 10, and another one 
(Fig. 8b) where, in addition, rj is reduced to 0.7, which has the 
effect of further increasing x a /x c (to ~ 1 x 10 3 ). We now sum- 
marize the distinguishing properties of these flow regimes. 

Outer region (x > x a ): ideal-MHD infall. This region ex- 
tends from x max = 10 2 , where the initial conditions are ap- 
plied, to x = x a , the location of the resolved ambipolar diffu- 
sion shock. This region does not differ substantially from the 
IMHD case, with the effects of ambipolar diffusion remaining 
minimal. The approximation b z k> b z> i ov/ (eq. [49]) is excellent, 
and it is also quite accurate to set [i = ajb z ^ const. The outer 
asymptotic solution remains roughly applicable for x > 10. Al- 
though the centrifugal force is still dynamically unimportant, 
there is some magnetic braking, especially forx < 1, and so the 
ratio j/m is not exactly constant. Over most of this region the 
infall speed is governed by self-gravity (m 3> m c ), but closer 
to the inner edge (x <~ 2), infall starts to be dominated by the 
central mass. 

The ambipolar diffusion shock (x w x a ). This shock marks 
the end of the IMHD regime. It can be resolved as a continu- 
ous transition, although it may contain a viscous subshock. Its 
most notable feature is the rapid increase in b z , which grows 
from essentially b z .i ow at the outer edge of the shock transition 
to ^z.high (eq- [50]) at the inner edge. Although the inequal- 
ity h <C x continues to hold, the db z /dx terms in the structure 
equations are important in this range. [It is worth noting in 
this connection that the sharp spike exhibited by the disk half- 
thickness curve at the location of the shock is a consequence of 
the large £> z -gradient term in eq. (31). The extreme narrowness 
of the spike (Ax <C \u\h) indicates, however, that the assump- 
tion of vertical hydrostatic equilibrium is not justified at that 
location. In reality, the enhanced magnetic squeezing at the 
shock will reduce the disk thickness to only a fraction of the 
amount indicated by the equilibrium curve over the fluid tran- 
sit time through the shock. This apparent overshoot is, how- 
ever, a highly localized phenomenon: the solution obtained by 
omitting the db z /dx term in eq. (31) has a very similar global 
structure.] Immediately behind the shock the inward accelera- 
tion is temporarily reversed, resulting in a thin overshoot layer 
where the flow undergoes a weak outward acceleration. The 
shock structure is basically the same as that found in CCK for 



a nonrotational flow: for our fiducial parameters, rotation is not 
yet dynamically significant at this location. 

To obtain an estimate of the value of x a (which can be used 
as an initial guess in the numerical solution), we make use of 
the fact that, in the infall region that is located just downstream 
from x a , b ns w b z w b z .hi g h- These approximate equalities im- 
ply x 3 w (rjtp 2 /m)(h/ 'a) 1 ! 1 . We evaluate this expression at x a 
by using the fact that, for x > x a , ideal MHD is approximately 
valid and hence m/tp is not very different from its initial value 
Ho, and that at this location m is already close to m c so that 
tjj w m c /fio = const. For the solutions shown in Figures 7 and 
8, the vertical squeezing at x a is controlled by magnetic pinch- 
ing, so h ~ 2x 4 a/ip 2 (see eq. [32]). These approximations yield 
an estimate of the shock location, 

x a w V2r]/fi Q , (58) 

which depends only on prespecified flow parameters. Analo- 
gous (albeit less simple) algebraic expressions can be written 
down for the cases in which the vertical squeezing is dominated 
by the central tidal force or by self-gravity. 

Since the magnetic flux contained within x a is essentially the 
same as the flux that initially threaded the enclosed mass, it 
also nearly equals the flux trapped in the central split monopole 
of the IMHD solution. Ambipolar diffusion does not destroy 
the accumulated flux but rather redistributes it between the ori- 
gin and the shock location. As the value of x a generally in- 
creases with the diffusivity parameter 77, the approximate flux 
conservation within x a implies that the vertical field component 
(~ ^z.high) behind the ambipolar-diffusion shock decreases with 
increasing r\ (in particular, it scales as r?~ 2 when x a is given by 
eq.[58]). 

A comparison of Figures 8a and 8b reveals that the decel- 
eration of the flow in the AD shock is stronger the lower the 
value of rj. This can be attributed to the inferred dependence of 
^z.high on 77, which results in higher magnetic tension and pres- 
sure forces downstream of the shock when r\ is reduced. If r\ 
becomes smaller than ^0.2 (with the other parameters remain- 
ing unchanged), then the localized decrease of the infall speed 
inward of x a becomes large enough to reduce w = x — u below 
1 . But w = 1 is the singular line of the AD system of equations, 
corresponding to the critical (sonic) speed (see Appendix B). A 
crossing of this curve signals a discontinuity in the flow param- 
eters, which in this case implies the presence of an unresolved 
viscous subshock at the downstream end of the resolved MHD 
shock transition (see Li 1998). 

AD-dominated infall (x c < x < x a ). Over most of the region 
between the two shocks the centrifugal force remains dynami- 
cally unimportant, so matter moves almost in free fall, with the 
gravitational force dominated by the central mass. The ions and 
field are only weakly coupled to the neutrals, with the ion radial 
speed being effectively zero (except in the vicinity of x c ). The 
vertical field component continues to be given by b z w b z ,high 
to a very good approximation. In many respects this zone re- 
sembles the innermost region of the nonrotating flow studied by 
CCK, especially when it is fairly extended (x a 3> x c ), as it is in 
the two solutions presented in Figure 8. In particular, b z oc x~ l , 
ip?zx 2 b z , and (over a somewhat narrower range) a oc w ocx -1 / 2 . 
If this region is wide enough, magnetic braking can reduce j to 
an essentially constant value, j w j v \. The mass is also well ap- 
proximated by its plateau value (see eq. [35]), m w m v \ rts m c . 



COLLAPSE OF ROTATING MAGNETIC CLOUD CORES 



13 




1e-1 



1e-3 1e-2 1.31e-2 1e-1 1e0 1e1 1e2 
x = r/Ct 

Expanded 
Scale 



1e1 



1e0 



Centrifugal 
Shock Region 



_ _ — * \ 



AD 
Shock 

I 



(b) 



b x 

z 



1 x . b x=\i//x 

■ \ r,s T 



-b. /b 

<t>,s z 



I3c 



u D /u 



1e-1 

1e-3 1e-2 1.31e-2 1e-1 1e0 1e1 1e2 
x = r/Ct 

FIG. 7. — Fiducial solution for an ambipolar-diffusion rotational collapse. 
The variations of the normalized radial infall speed —u, surface density <r, to- 
tal mass m, mass accretion rate m, specific angular momentum j, and the disk 
half-thickness h as functions of the similarity variable x are plotted in Fig. 7a, 
and the corresponding variations of the z component (b z ) and the surface ra- 
dial (b r , s = 4>/x 2 ) and azimuthal (fo^,,) components of the magnetic field, as 
well as of the normalized radial component of the ion-neutral drift velocity 
(ud = Vo.r/C), are shown in Fig. 7b. The diffusivity parameter is t] = 1. The 
other parameter values are vo = 0.73, A = 3, uq = — 1, fj,Q = 2.9, a = 0.08, and 
<5=1. The ambipolar-diffusion and centrifugal shocks are located at x a = 0.41 
and x c = 1.3 X 10~ 2 , respectively. The central mass is m c = 4.7, and the 
plateau values in the AD-dominated infall region are m p \ m 4.9 and j p \ 0.26. 
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FIG. 8. — (a) Same parameter values as in Fig. 7, except that vo is re- 
duced from 0.73 to 0.18. In this case x a = 0.46, x c = 1.7 X 10~ 3 , m c = 6.0, 
m p i = 6.1, and j p \ f» 0.1. (b) Same as Fig. 8a, except that the diffusivity pa- 
rameter r\ is reduced from 1 to 0.7. In this case x a = 0.33, x c = 2.6 X 10 — 4 , 
m c = 6.00, m p i = 6.03, and j p \ « 0.04. The horizontal scale in the vicinity 
of the centrifugal shock has been expanded to show details of the postshock 
transition zone. 
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As we point out in § 4, m c will be reduced below the value of 
m upstream of the centrifugal shock if (as expected) the rota- 
tionally supported disk drives a wind from its surfaces. The nu- 
merical value of the reduction factor depends, however, on yet 
another model parameter (see eq. [C6]). Therefore, to simplify 
the discussion, we have chosen not to incorporate the effect of 
the wind mass loss into the solutions that we present. 

The centrifugal force starts to become important near the 
downstream end of this zone: after it comes to exceed the grav- 
itational force, it triggers the formation of the centrifugal shock 
at the point x c w j v \ jm c . In contrast to the IMHD solution pre- 
sented in § 3.2.1, one cannot approximate j v \ by its initial value 
in this case: the braking action induced by the magnetic field 
amplification in the AD shock completely invalidates this ap- 
proximation. A comparison of Figures 8a and 8b illustrates this 
point. When r\ is reduced from 1 (Fig. 8a) to 0.7 (Fig. 8b), the 
location of the AD shock changes by a factor of ~ 0.7 and the 
value of the central mass remains essentially the same: this is as 
expected from equations (58) and (35), respectively. However, 
the decrease in rj has a very pronounced effect on the location 
of the centrifugal shock: as r\ is reduced by a factor of <~ 1 .4, x c 
decreases by a factor of ~ 6.7. As we now demonstrate, this can 
be traced to the sensitive dependence of the magnetic braking 
term on rj, which leads to a strong reduction in j downstream 
from the AD shock for even a moderate decrease in rj. 

We first estimate bfo by substituting b z from equation (50) 
into equation (24) and using equation (25), 

b^ « ~2a^jx-\\ +2aw)- 1 . (59) 
Equation (23) then yields 

^ « -(x 2 /m)b z b 4> , s (60) 

^2ar]- 1 j(a/h) 1/2 (l+2aw)- 1 . (61) 

Using again the approximations b v w b z and h rts 2a/b 2 s that 
were employed in deriving the above estimate of x a , equation 
(61) becomes 

dj am j 

— « -5 — j(l +2aw) . (62) 

ax r) L x 

Using the approximations w w (2m/x)'/ 2 and m w m c , and as- 
suming that 2aw 3> 1, we get 

d -L*(rnM /2 ri- 2 x- l/2 j. (63) 

This equation can be integrated to give 

7 = 7iexp{(m c /2) 1 /VV/ 2 } , (64) 

where j\ is a constant. This result is valid between x a and x c and 
can be used to estimate the latter. For x c <C x a , the value of j 
at x c may be identified with j v \ and can be approximated by j\ . 
Assuming also that no significant magnetic braking takes place 
for x > x a (consistent with our analysis of the IMHD solution 
in § 3.2.1), we can write j v \ w (vo/A)m p i w (vo/A)m c . Using 
equation (58) to evaluate x a , we finally obtain 
2 

x c «^|m c exp{-(2 3 / 2 r,- 3 m c / M0 ) 1/2 } . (65) 

We have verified that this estimate (which depends on the ratio 
x a /x c being 3> 1) yields the value of x c to within 30% for both 
of the solutions shown in Figure 8. The strong dependence of 
x c on 77 that is exhibited by these solutions can be understood 
from the form of the expression (65): it indicates that the ef- 
fect of even a small variation in 77 is magnified on account of 



its presence inside an exponential function. Although the pre- 
cise form of this expression may be expected to differ in cases 
where some of the approximations we utilized do not apply, the 
sensitive dependence of x c on rj is probably a general feature of 
the model. 

The centrifugal shock (x w x c ). The shock represents a dis- 
continuity in the variables u, a, and b z through which the flow 
changes from the nearly undiluted free fall to subsonic accre- 
tion. We approximate the shock as a discontinuity in the vari- 
ables u and a, which is governed by the isothermal-shock jump 
conditions (Appendix B . 1 ), followed by a very thin (but numer- 
ically resolvable) layer where b z increases. In this picture, the 
ions are decoupled from the neutrals inside the subshock, so 
their radial speed in the frame of the shock, and hence (by mag- 
netic flux conservation, see eq. [A18]) B z , remain unchanged 
across the subshock even as \V r \ and S vary. However, the in- 
crease in the drift speed \Vo. r \ across the subshock enhances the 
collisional drag force immediately behind it and causes a local 
recoupling of ions and neutrals; this, in turn, leads to a rapid de- 
cline in I Vj >| and a consequent increase in B z . Just inward of the 
shock there is a narrow transition zone where the flow settles 
down to its asymptotic Keplerian structure. Within this layer 
the variables u, a, and b z adjust rapidly (with some oscillations 
and overshoots) to their asymptotic values. This implies, in par- 
ticular, that, just as in the case of the AD shock, the derivative 
terms hdb z /dx in the constituent equations cannot be neglected 
in this zone. The surface density in this region is significantly 
larger than a in its immediate vicinity, so the layer represents 
a massive ring lying just outside the main body of the disk (re- 
sembling the situation in the fiducial IMHD solution considered 
in § 3.2.1). In the solution shown in Figure 7, the width of the 
ring is ~ 10% of x c and it encompasses ~ 8% of the total mass 
of the disk within x c ; the relative size of the ring in the solutions 
displayed in Figure 8 is, however, smaller. The variables ip and 
j do not change appreciably across this region. 

The Keplerian disk (x < x c ). Inward of the postshock tran- 
sition layer, the flow approaches the form of a magnetized Ke- 
plerian accretion disk as described by equations (51)— (57). For 
the solutions shown in Figures 7, 8a, and 8b, the total mass 
of the rotationally supported disk (including the ring at x c ) is, 
respectively, ~ 5%, <~ 2%, and <~ 0.4% of m c . 

3.3.2. Fast Rotation 

As in the corresponding IMHD case, a high value of vo and 
comparatively low values of a and S result in the centrifugal 
shock being located within the self-gravity-dominated region 
(where m ^> m c ). The value of x c is large enough for ideal MHD 
to still be applicable (so tp w xwb z ). Since the inflow speed is 
strongly reduced in the dense disk that forms behind the cen- 
trifugal shock, the ambipolar-diffusion time becomes shorter 
than the accretion time and an AD transition front is established 
not far from x c . In fact, as — u is small enough for the flow to 
remain subsonic (x — u < 1) in this region, the IMHD-to-AD 
transition is gradual rather than sharp as in the other cases con- 
sidered in this subsection, where it takes the form of an AD 
shock. Since the centrifugal shock occurs within the ideal- 
MHD region, we impose the same jump conditions at x c as in 
the corresponding IMHD case (see § B.2). 

A representative solution is shown in Figure 9. It exhibits the 
same general features (a non-Keplerian outer region, a small 
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FIG. 9. — Fast-rotation solution for an ambipolar-diffusion rotational col- 
lapse. The variations of the normalized radial infall speed — u, surface density 
<t, total mass m, specific angular momentum j, and z component of the mag- 
netic field b, are plotted as functions of the similarity variable x. The model 
parameter values are the same as in the fiducial case (Fig. 7), except that the 
azimuthal- velocity parameter vo is increased from 0.73 to 1.5 and the braking 
parameter a is increased from 0.08 to 0. 1 (so they match the corresponding 
parameters in the ideal-MHD solution depicted in Fig. 4). The derived values 
of the centrifugal-shock radius (x c = 1.5) and of the central mass (m c = 0.5) 
are similar to those of the corresponding ideal-MHD solution. 

central mass, a backflowing layer behind the centrifugal shock) 
as the fast-rotation IMHD solution (see Figs. 4 and 5). These 
characteristics are basically a consequence of the fast initial ro- 
tation, and the strong similarity with the corresponding IMHD 
case shown in Figure 4 can be understood from the fact that the 
centrifugal shock occurs outside the AD regime. 

3.3.3. Strong Braking 

In § 3.2.3 we have found that ideal-MHD flows with strong 
braking can lose all their angular momentum at a finite distance 
from the origin. Here we consider whether a similar situation 
can occur in the presence of ambipolar diffusion. 

When 7 = there are two inner asymptotic solutions that are 
mathematically self-consistent and that can be matched to the 
outer asymptotic solution by integrating the system of differ- 
ential equations. Correspondingly, two kinds of solutions are 
found: one that is almost IMHD in nature, and another that is 
AD dominated. 

The IMHD-like solution resembles the one obtained in 
§ 3.2.3 in that b z is everywhere well approximated by the root 
b Z: \ \N of equation (48), there is no AD shock, and the inner 
asymptotic solution is given by equations (43)-(47) (but with 
\i possibly differing from ^ on account of the small, but fi- 
nite, ambipolar diffusion). Outside the point xj the azimuthal 
field component can differ significantly from its IMHD value, 
but once complete braking is accomplished and both v and Z?^ , 
effectively vanish, this solution basically describes an IMHD 
nonrotating inflow and, in particular, predicts the formation of 
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FIG. 10. — Strong-braking solution for an ambipolar-diffusion rotational 
collapse. The variations of the normalized radial infall speed — u, surface den- 
sity <x, total mass m, specific angular momentum j, and z, component of the 
magnetic field are plotted as functions of the similarity variable x. The model 
parameters are rj = 0.5, v = 1, a = 10, and <5 = 10, with the remainder having 
the same values as in the fiducial case (Fig. 7). Note that j attains a tiny, but 
finite, plateau value in the innermost (free-fall) region (j p \ = 2 X 10 -5 ). The 
AD shock is located at x a = 0.24, and the central mass is m c = 5.9. 



a central split-monopole field. As we remarked at the begin- 
ning of § 3.3, such a configuration is unlikely to occur in real 
systems, so we henceforth concentrate on the more plausible 
alternative case. 

In the diffusive solution, the flow passes through an AD 
shock at a point x a (which can be estimated by eq. [58]) that 
typically exceeds xj. Inward of the ambipolar-diffusion shock, 
the poloidal field components have the same form as in the fidu- 
cial solutions, b rs w b, k> £> z .high- As in the corresponding IMHD 
case, the inner asymptotic solution satisfies j = b^^ w and de- 
scribes a free fall onto the central mass, which, however, in this 
case is not magnetically diluted but rather resembles the situa- 
tion in the AD-dominated infall region of the fiducial solutions. 
It is given by 



m ■ 



■ m - 
-u - 
a ■■ 



b r . s 



(2m c /x) 1 / 2 , 
(m c /2x) 1/2 , 
b z = l-^m^rj-^h-^x-' 

h = hx 3/2 . 



(66) 
(67) 
(68) 
(69) 
(70) 
(71) 



The power-law scalings of the flow variables with x are the 
same as in the asymptotic solution obtained by CCK for the 
nonrotating AD collapse. However, because of the magnetic 
pinching effect included in the above solution, the expression 
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for the coefficient h is more involved: 



h = 



(12D/VG-G) l / 2 -VG 
2^6 

D=N l/4 /rj, 
N = 2/m c , 

G=(2F) l ' 3 -W{3/F) 1 ' 3 , 

F = 9D 2 + [3(27D 4 + 256JV 3 )] 1/2 



(72) 

(73) 
(74) 
(75) 
(76) 



with the expression for h representing the real, positive solution 
of the quartic equation h 2 +Dh x / 2 —N = 0. 

An example of a strongly braked solution is shown in Figure 
10. A curious feature of the flow is the value of j as x — > 0, 
which assumes a nonzero (albeit rather small) plateau value 
Opi ~ 10~ 5 ). This puzzling behavior can be attributed to the 
very small magnitude of ip in this limit, which causes \b^J[ (eq. 
[24]) and therefore \dj/dx\ (eq. [23]) to be small, more than 
compensating for the increase in b z behind the AD shock. The 
low value of the angular-momentum derivative in turn accounts 
for the asymptotic near-flatness of the j curve. (In contrast, 
the large central value of ip associated with the split-monopole 
field in the near-IMHD solution guarantees that, in this case, 
\b,p :S \ and \dj/dx\ remain large enough near xj to make j ac- 
tually vanish there.) The value of j p i is highly sensitive to the 
magnitude of the parameter r\\ it decreases by <~ 10 orders of 
magnitude when t] is decreased from 0.5 to 0.3. (The strong 
sensitivity of x c to the value of this parameter has already been 
noted in connection with the fiducial solutions; see eq. [65].) 
However, for practical purposes, j p i can be regarded as being 
effectively zero. In particular, the implied location of the cen- 
trifugal radius, x c w j v 2 /m c , corresponds to such a small radius 
for any meaningful value of t that it would invariably lie well 
inside the central protostar. 

4. DISCUSSION 

The model presented in this paper provides a physical frame- 
work for a detailed study of the development of star-disk sys- 
tems in the collapse of rotating, magnetic molecular-cloud 
cores. Perhaps the most noteworthy feature of the solutions 
presented in § 3 is that, for representative values of the initial 
cloud rotation and magnetization and with plausible estimates 
of the magnetic field diffusivity and of the strength of magnetic 
braking, they predict the formation of circumstellar disks whose 
properties compare quite well with those of real protostellar 
systems. In particular, the disks obtained for our chosen fidu- 
cial parameters are Keplerian, having a mass that is only a small 
fraction (< 10%) of the central mass, as has been inferred for 
most protostellar disks in nearby dark clouds (see § 2. 1). For the 
representative AD solution shown in Figure 7 (corresponding to 
r\ = 1, Vo = 0.73, a = 0.08, and S = 1), the disk mass Md is less 
than 5% of the central mass M c . For comparison, if the model 
parameters are slightly changed (to r\ = 1, Vo = 0.77, a = 0.1 1, 
and 5 = 0.7), one finds Md /M c w 9%. The most reliable proto- 
stellar disk velocity curves have been obtained in Class II YSOs 
(i.e., T Tauri stars, of typical ages <~ 10 6 yr), which, in contrast 
to younger protostars, are no longer strongly obscured by dusty 
envelopes (identified with the contracting cloud cores). These 
are also the systems with the most accurate determinations of 
the disk mass (e.g., Mundy et al. 2000): the massive envelopes 
surrounding Class YSOs prevent their disk masses from being 
pinned down, and even in the case of Class I sources the values 



of Md are not as well determined (although they appear to be 
similar to the disk masses in Class II systems). It has occasion- 
ally been suggested in the literature that protostellar disks may 
not be Keplerian during the early (in particular, Class-0) phase 
of their evolution, when the ratio of the central mass to that of 
the surrounding envelope is still small. Our results demonstrate 
that this need not be the case. In fact, in the self-similar solu- 
tions that we derive, if the mass of the rotationally supported 
disk is small in comparison with the central mass at any given 
time, it will be so at all times. 

The observed rotationally supported disks typically have 
sizes < 10 2 AU (see § 1), which are reproduced by our typi- 
cal solutions. For example, the derived disk size (given by the 
value of the centrifugal radius r c ) 10 s yr after the start of the 
collapse (which corresponds to the Class-I evolutionary phase) 
is ~ 52 AU for the fiducial solution presented in Figure 7 and 
<~ 130 AU for the other AD solution mentioned in the preceding 
paragraph. There have been reports in the literature of T Tauri 
circumstellar disks that exhibit Keplerian rotation on somewhat 
larger scales [e.g., GM Aur, with r c > 200 AU (Dutrey et al. 
1998), and DM Tau, with r c > 600 AU (Guilloteau & Dutrey 
1998). It is quite possible that these disks have achieved their 
large sizes already during their early evolutionary stages. For 
example, they may have possessed a comparatively high value 
of the diffusivity parameter 77, to which the predicted magni- 
tude of r c is particularly sensitive (see eq. [65]). Alternatively, 
they may have formed from collapsing cores with compara- 
tively large initial rotations and relatively inefficient braking 
(see § 3.3.2). It is worth noting, though, that in the self-similar 
model, disks with larger values of x c tend to have higher disk- 
to-star mass ratios and therefore exhibit progressively greater 
departures from a Keplerian rotation law. It is, however, also 
conceivable that at least part of the cause for large observed disk 
sizes in these highly evolved (and no longer strongly accreting) 
systems is post-formation viscous spreading of the disk (e.g., 
Lin & Pringle 1990; Hartmann et al. 1998). It is interesting to 
note in this connection that the radial scaling of the disk sur- 
face density predicted by our self-similar solution (S oc r~ 3 / 2 
for small r; see eq. [54]) is the same as that inferred for the 
"minimum mass" solar nebula (Weidenschilling 1977; Hayashi, 
Nakazawa, & Nakagawa 1985). Although this agreement is en- 
couraging, it is not unique to this model: in particular, the same 
scaling is predicted for a self-similar Keplerian disk with "a 
viscosity" (Tsuribe 1999) and for a self-similar self-gravitating 
disk with a gravitational instability-induced effective viscosity 
(Lin & Pringle 1987). 

Our model disks exhibit a surface-density enhancement near 
their outer edges (see § 3.2.1 and § 3.3.1): for the fiducial AD 
solution shown in Figure 7, this ring has a mass <~ 0.1 Md and 
a width ~ 0.13 r c . We are not aware of any observational evi- 
dence as yet for the existence of such rings in protostellar sys- 
tems, but there have been several recent reports of well-defined 
dust rings with similar radii (~ 10 2 AU) and fractional widths 
around older, Vega-like stars (e.g., Schneider et al. 1999; Dent 
et al. 2000; Koerner, Sargent, & Ostroff 2001). Although the 
apparent confinement of the latter rings is usually attributed to 
the presence of solid "satellites" (e.g., planets), it is interesting 
to speculate that their origin could be traced to the E enhance- 
ments predicted by our solutions. 

An important implication of the result that typical rotating 
inflows tend to produce a central mass that is surrounded by 
a Keplerian disk is that the angular momentum problem (see 
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§ 1) for the modeled YSO is basically resolved. In particular, 
angular momentum transport can be sufficiently efficient to al- 
low most of the inflowing mass to end up (with effectively no 
angular momentum) at the center, with the central mass dom- 
inating the dynamics well beyond the outer edge of the disk 
even as the inflow is still in progress. Our solutions reveal that 
the AD shock, even though it is usually located well outside the 
region where the centrifugal force becomes important, helps to 
enhance the efficiency of angular momentum transport through 
the magnetic field amplification that it induces. In practice, the 
ultimate value of the protostellar angular momentum would be 
determined by the braking mechanism that enables the YSO 
to reduce its angular velocity below the centrifugal-equilibrium 
value that characterizes the accreted disk matter (e.g., Mestel & 
Spruit 1987; Konigl 1991; Popham 1996). 

The formation of Keplerian circumstellar disks in our "stan- 
dard" solutions is a consequence of efficient magnetic brak- 
ing, which allows the mass injected at large distances to 
reach the center rather than pile up in the disk (as it does in 
the fast-rotation solutions discussed in § 3.2.2 and § 3.3.2). 
The disk in these cases therefore has little influence on 
the mass accretion rate onto the central YSO, M c = 1 .59 x 
10 _6 (r/10K) 3 / 2 m c M Q yr _1 , which is determined on scales 
much larger than r c . In fact, as we discuss in § 3.1, the eigen- 
value m c is typically also insensitive to the magnetic structure 
of the flow and is often well approximated by the "plateau" 
value OTpi (eq. [35]) that depends only on the initial condi- 
tions. (For our fiducial parameters, m v \ = 6, implying M c w 
9.5 x 10~ 6 M Q yr _1 .) The likely production of a centrifugally 
driven wind that transports mass away from the disk surfaces 
(see discussion below) will, however, act to reduce the value of 
m c : as we estimate in Appendix C, the reduction factor could 
be as large as ~ 3. 

The presence of rotation also does not strongly affect the 
"revitalization" of ambipolar diffusion behind the AD shock, 
which, as we noted in § 1, can go a long way toward resolv- 
ing the magnetic flux problem. This conclusion applies even in 
the fast-rotation case (see § 3.3.2), where the AD shock is re- 
placed by a more gradual transition occurring within the disk. 
The asymptotic x — > similarity solution for the disk implies 
that, at any given time, the magnetic flux scales as r 5 / 4 (see 
eq. [56]), which represents a slower decrease with r than in the 
nonrotating case (in which $ oc r; see eq. [70]). However, in 
reality, the amount of magnetic flux that ends up threading the 
protostar would depend on the detailed flux transport and de- 
struction mechanisms that operate in the inner disk and the star, 
which may involve Ohmic resistivity (e.g., Li & McKee 1996), 
ambipolar diffusion (e.g., Desch & Mouschovias 2001), and re- 
connection (e.g., Goodson & Winglee 1999). 

The quasi-equilibrium disk configurations described by our 
solutions make it possible to investigate the stability properties 
of the modeled circumstellar disks. Although a full pursuit of 
this topic is beyond the scope of this paper, we comment briefly 
on the basic MHD and gravitational instabilities that may affect 
the disks. The most relevant MHD instability is evidently the 
magnetorotational one, which is expected to develop rapidly in 
any differentially rotating disk that is threaded by a subther- 
mal magnetic field (see Balbus & Hawley 1998 for a review). 
This instability has been shown to generate turbulence that can 
remove angular momentum from the inflowing gas, so the de- 
termination of whether the disk is susceptible to this instability 
also has implications to the self-consistency of our identifica- 



tion of magnetic braking as the dominant angular-momentum 
transport mechanism in the disk. 

As it turns out, our diffusive Keplerian disk models are by 
and large magnetorotationally stable. Physically, this is be- 
cause the coupling between the matter and the field is gener- 
ally too weak to allow the instability to grow. Specifically, a 
necessary condition for this instability to operate in a weakly 
ionized, Keplerian disk is r n ; < fij^ 1 , where £Ik is the Keplerian 
angular velocity (Blaes & Balbus 1994). Using our self-similar 
variables, this criterion can be written as (ax 3 /m c hrj 2 Y I 1 > 1. 
Since the Keplerian disk region in all of our typical solutions 
nearly coincides with the region where the inner asymptotic so- 
lution (eqs. [5 1]— [57]) is applicable, we can use the latter to 
evaluate the left-hand side of the above inequality; it is found 
to equal 2/35. [The asymptotic solution satisfies the weak-field 
condition V^./C 2 < 1 postulated in the derivation of the insta- 
bility criterion; indeed, the left-hand side of the latter inequality 
(which equals b 2 h/a) vanishes (ocx 1 / 2 ) as x — > 0.] But the pa- 
rameter S (defined in § 2. 1) is expected to be > 1, and therefore 
the instability condition T n iQ K < 1 is not satisfied. Interestingly, 
our fiducial IMHD solution is also magnetorotationally stable. 
In this case the strong magnetic squeezing prevents the unsta- 
ble perturbations from fitting inside the disk. Formally, one re- 
quires the critical wavelength A c w 3.6(Va z /^k) to be less than 
the disk thickness 2H for the disk to be unstable. Using the in- 
ner asymptotic solution given by equations (37)-(42), with the 
help of which we first verify that V^/C 2 again tends to zero 
as x — > 0, this condition translates into 1.8(m c /2/x 4 x)'/ 2 < 1. 
For our typical parameters this requires x > 0. 14, but this lower 
bound exceeds x c (the approximate disk outer boundary) for this 
solution. 

The above results support the choice of magnetic braking 
as the main angular momentum transport mechanism in our 
disk model. This mechanism remains effective even in the 
presence of ambipolar diffusion because the neutral-ion col- 
lision time T n ;, while being comparatively long, is shorter than 
the magnetic braking time Tt, r . In a Keplerian accretion disk 
Tb r w r/\V r \ and V r is typically much smaller than the Keple- 
rian speed; therefore T^jr^ can be <C 1 even if t^^Ik exceeds 
1 . This conclusion continues to hold when the angular momen- 
tum transport in the disk (and hence the accretion time r/|V r |) 
is controlled by a centrifugally driven wind rather than by the 
propagation of torsional Alfven waves. As we show in Ap- 
pendix C, the basic formalism of our model is unchanged in 
this case, and the parameter 5 (which is now determined di- 
rectly from the centrifugal wind solution) is again typically 
small enough for the linear stability criterion t^SIk > 1 to be 
satisfied over most of the disk column density. The argument 
against instability is reinforced by indications from numerical 
simulations (e.g., Hawley & Stone 1998) that MRI-driven tur- 
bulence and the resulting angular-momentum transport are sig- 
nificant only when T^flx -C 1. We note in this connection that 
even a disk in which the bulk of the material is weakly coupled 
is expected to possess well coupled (Tm^K < 1) surface layers 
from which the centrifugal wind is launched (see Appendix C). 
It is, however, likely that the critical wavelength for instability 
exceeds the density scale height in the strongly coupled regions 
(see Wardle & Konigl 1993 and Li 1996), so the surface layers 
should typically also be magnetorotationally stable. 

The basic gravitational instability that we consider involves 
fragmentation induced by the disk self-gravity. The rele- 
vant stability condition against axisymmetric perturbations is 
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given by the Toomre criterion (2Toomre = Ck/tiGE > 1 (Toomre 
1964), where k = [r~ 3 (dJ 2 /dr)] 1 / 2 is the epicyclic frequency. 
Although this criterion is modified somewhat by magnetic ef- 
fects (Shu & Li 1997; R. Krasnopolsky & A. Konigl, in prepa- 
ration), the corrections remain small for nearly Keplerian disks, 
especially in the AD regime, and we therefore neglect them 
here. We have verified that the Toomre criterion is well satis- 
fied for most of our rotationally supported disk solutions, the 
only exceptions being the outermost (backfiow) regions of the 
fast-rotation configurations (Figs. 4, 5, and 9), where the effect 
of self-gravity is comparatively large. This result is significant 
in that it suggests that in typical cases protostellar disks may 
not fragment even during the early phases of their evolution, 
when the central mass is still small. In our model this behavior 
can be attributed to the fact that, for representative parameters, 
the central gravity already dominates at the location of the cen- 
trifugal shock, so that the disks are nearly Keplerian from the 
start. Inasmuch as the inequality ^Toomre < 1 signals the on- 
set of self-gravity-induced angular-momentum transport in the 
disk (e.g., Lin & Pringle 1987), this result also serves to jus- 
tify the neglect of such a transport mechanism in our model. 
We note, however, that the collapsing core could in principle 
be unstable to fragmentation before point-mass formation, or, 
following PMF, in the region outside the rotationally supported 
disk (e.g., Boss 2000). It is also conceivable that fragmenta- 
tion might be initiated even before a dynamically unstable core 
is formed (possibly triggered by an AD-mediated instability in 
the magnetically supported parent cloud; see Zweibel 1998 and 
Indebetouw & Zweibel 2000). 

Another potential angular-momentum transport mechanism 
that was not included in our model involves centrifugally driven 
disk outflows. For a cold, thin, Keplerian disk that is threaded 
by open magnetic field lines, such an outflow could be launched 
if the angle between the meridional projection of the magnetic 
field and the rotation axis exceeds 30° (Blandford & Payne 
1982). This condition is satisfied by our asymptotic disk solu- 
tions — in particular, the AD inner asymptotic solution implies 
a radially constant inclination angle <~ 53° (see eq. [56]) — so 
our model disks can be expected to give rise to outflows of this 
type. This result is encouraging in view of the fact that cen- 
trifugally driven disk outflows are a leading candidate for the 
origin of the bipolar outflows that are frequently observed to 
emanate from YSOs (see, e.g., Konigl & Pudritz 2000 for a 
review). A centrifugally driven disk outflow can be incorpo- 
rated into our model by taking account of the fact that the wind 
would be quasi-steady on the accretion time scale and should 
thus be well approximated by a steady-state formulation. The 
condition that the centrifugal outflow control the vertical angu- 
lar momentum transport from the disk is equivalent to the re- 
quirement that the steady-state wind solution pass through the 
Alfven critical point. This requirement, in turn, serves to fix 
the value of B^, s (or, equivalently, of the parameter 5) in the 
disk solution. Besides angular momentum, the wind also car- 
ries away mass, resulting in a progressive decline in the mass 
accretion rate with decreasing disk radius and leading to lower 
values of m and m at the center. Interestingly, the radial scaling 
of the disk magnetic field (oc r -5 / 4 ; see eq. [55]) implied by 
the asymptotic AD disk solution for any given value of t is the 
same as in the Blandford & Payne (1982) radially self-similar, 
steady-state wind solution. This suggests that the latter solution 
can be used in conjunction with our disk model to fix the val- 
ues of 8 and of the m c reduction factor in the asymptotic disk 



solution. We discuss the implementation of this scheme in Ap- 
pendix C. Inasmuch as the incorporation of a disk wind does 
not qualitatively change our basic results, we have not included 
its effects in the model described in § 2 so as not to unduly 
complicate the presentation. It is nevertheless important to rec- 
ognize that our solutions indicate that the collapse of rotating, 
magnetic cloud cores naturally gives rise to disks with mag- 
netic field configurations that are conducive to the launching of 
centrifugally driven winds. This conclusion is supported by the 
MHD core-collapse numerical simulations of Tomisaka (1998, 
2000, 2002), in which disk outflows that transport angular mo- 
mentum through material and magnetic stresses are observed 
to form. (The possibility that a disk formed in the collapse 
of a magnetized rotating cloud core could give rise to a cen- 
trifugally driven wind was recognized also by Contopoulos & 
Sauty (2001). However, in contrast with our model, in their pic- 
ture the angular momentum transport in the disk is dominated 
by turbulent viscosity and the outflow is concentrated near the 
outer disk radius; furthermore, they infer that the mass outflow 
rate depends sensitively on the efficiency of magnetic braking 
in the precollapse core.) 

We have also presented solutions (§ 3.2.3 and § 3.3.3) in 
which the magnetic braking is sufficiently efficient to prevent 
the formation of a rotationally supported disk from the collapse 
of a rotating cloud core. Although we argued that the param- 
eter combinations characterizing these solutions are not typi- 
cal, it is conceivable that such situations could be realized. A 
possible indication that this may sometimes be the case is pro- 
vided by the detection of slowly rotating YSOs that show no 
evidence (from excess near- and mid-IR emission) for the pres- 
ence of a circumstellar disk (Stassun et al. 1999, 2001). These 
observations are puzzling if one considers disk accretion to be 
a ubiquitous feature of star formation and interprets slow stel- 
lar rotation as due to a star-disk interaction (e.g., Konigl 1991; 
Popham 1996). The properties of these YSOs are, however, 
readily understood if they are the product of a strongly braked 
core collapse. 

5. SUMMARY 

We study the collapse of rotating, magnetic molecular-cloud 
cores with the help of a self-similar model. This model gener- 
alizes the work of CCK on nonrotating cloud cores by incorpo- 
rating the effects of rotation and magnetic braking (while also 
taking into account ambipolar diffusion and its dependence on 
the magnetic tension force as in the original model). We fo- 
cus on the evolution after point-mass formation, noting that the 
pre-PMF collapse has already been studied by numerical sim- 
ulations. Our approach is motivated by previous findings that 
simulation results for the pre-PMF collapse of rotating cores 
and for the pre- and post-PMF collapse of nonrotating cores 
are well approximated by self-similar solutions. Our semiana- 
lytic scheme allows us to examine the full range of expected 
behaviors and their dependence on the physical parameters. 
We present and analyze solutions for rotating but nonmagnetic 
collapsing cores, reproducing the results of Saigo & Hanawa 
(1998), as well as for magnetic cores in the ideal-MHD and 
ambipolar-diffusion limits. In these two regimes, we distin- 
guish between fiducial solutions, obtained for typical parameter 
values that correspond to moderately fast initial cloud rotation 
rate and moderately strong magnetic braking, and limiting cases 
of (1) fast rotation and (2) strong braking. Our results can be 
summarized as follows: 
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• For representative parameter values, we obtain solutions 
that describe rotationally supported circumstellar disks 
of masses and sizes that are consistent with observa- 
tions of YSOs. Our model thus makes it possible, for 
the first time, to study the formation of protostellar disks 
in the context of a realistic scenario of star formation in 
magnetically supported, weakly ionized molecular cloud 
cores. The outer boundary of the nearly stationary ac- 
cretion disk roughly coincides with the location of the 
centrifugal shock, which typically occurs inward of the 
ambipolar-diffusion shock where the magnetic field de- 
couples from the matter. 

• Our solutions indicate that it is quite possible that T Tauri 
(Class II) protostellar systems, whose disk masses are 
typically inferred to be < 10% of the central mass, have 
had a similarly low disk-to-star mass ratio also during 
their earlier (Class-0 and Class-I) evolutionary phases. 

• The disk configurations that we obtain exhibit surface- 
density enhancements near their outer boundaries (en- 
compassing < 10% of the disk mass and radius for typ- 
ical parameters). We speculate that these features may 
be the precursors of the rings detected in some Vega-like 
systems on scales < 10 2 AU. 

• Our models allow us to elucidate the interplay between 
the ambipolar diffusion and magnetic braking processes. 
In particular, the magnetic field enhancement behind the 
AD shock increases the efficiency of angular-momentum 
removal from the disk. In the limiting case of strong 
braking, essentially all the angular momentum is re- 
moved well before the inflowing gas reaches the center. 
Such systems may correspond to slowly rotating YSOs 
that show no evidence of a circumstellar disk. 



• The inferred magnetic-field structures imply that the 
disks could drive centrifugal outflows over much of their 
radial extent. We show that the steady-state, radially 
self-similar wind solution of Blandford & Payne (1982) 
can be naturally incorporated into the asymptotic AD 
disk solution, making it possible to study the effects 
of wind angular-momentum and mass removal from the 
disk and to better constrain the relevant parameters. 

• The derived nearly-Keplerian disk configurations appear 
to be immune to both the magnetorotational instability 
and to self-gravity-induced fragmentation. Only in the 
limiting case of rapid initial rotation (which yields com- 
paratively massive and largely non-Keplerian disk solu- 
tions) is the Toomre Q parameter < 1 anywhere in the 
disk. 

Despite the various simplifications that needed to be em- 
ployed to ensure self-similarity, our solutions likely capture the 
main qualitative aspects of the collapse of rotating cloud cores, 
and it is conceivable that our numerical estimates are also not 
too far off. Further progress in the study of protostellar disk 
formation could be achieved by means of numerical MHD sim- 
ulations. 
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APPENDIX 
A. DISK EQUATIONS 

Assuming axisymmetry and isothermality, we write down (using cylindrical coordinates) the mass, radial momentum, and angular 
momentum conservation relations as well as the vertical hydrostatic equilibrium condition: 

dp 1 8 8 

-£ + --(rpV r ) = --(pV z ) , (Al) 
at r or 8z 

dV r , „x/ _ r^OP , „ V , B z 8B r 8 ( B\ \ 1 d , D ^2 T , 8V r 



where g r and g z are, respectively, the radial and vertical components of the gravitational field. 

We now integrate these equations over z. We assume that the disk is geometrically thin [half-thickness H(r) <C r] and aim to get 
expressions that are valid to order (H/r) 2 . On account of the disk thinness, we approximate the density, radial velocity, azimuthal 
velocity, and radial gravity as being constant with height. We thus equate the column density S = pdz to 2Hp and drop the last 
term on the right-hand side in both equations (A2) and (A3). Equation (Al) becomes 

<9S + 1_<9 (r£V) _ l_8M« (A5) 

8t r 8r 2irr dr ' 

where the last term on the right-hand side represents the mass flux in a disk outflow [of total outflow rate M w (r) within a radius r]. 
We neglect this term in the calculations presented in the body of the paper, but we return in § C to consider the effect of having 
M w / o. 
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Using the solenoidal condition on the magnetic field, [dB z /dz — —r~ l (d/dr)(rB r ), which in a thin disk allows us to treat B z 
as being constant with height whenever it is not explicitly differentiated with respect to z] as well as the assumed field symmetry, 
equation (A2) integrates to 



yi B z B rs 



Or r 2tt 



2H 



d (B 



dr V 8tt 



(A6) 



Taking account of the adopted vertical mass distribution, we replace the integral in the last term of equation (A6) by an an integral 
over a finite interval, J^g) r \ dz: 



(A7) 



(see Lovelace, Romanova, & Newman 1994). We assume B r (r,z) — B ns (r) \z/H(r)\ , and similarly for B^. This choice is motivated 
by the field configuration derived in the case of a rotationally supported thin disk in which the field is comparatively well coupled 
to the matter (e.g., Wardle & Konigl 1993). Although the disks obtained in our AD collapse solutions evidently correspond to the 
weak-coupling case (see § 4 and Appendix C), for which this approximation is no longer adequate (e.g., Li 1996; Wardle 1997), we 
nevertheless adopt the above scalings for the sake of definiteness. We emphasize, however, that none of the dominant terms in the 
equations that we solve depends on the details of the vertical variation of these field components. We thus get 



'dH\ 

~d7) 



(A8) 



Approximating g r by GM(r)/r 2 [with M(r) w M c if the central mass dominates] and substituting J = rV^, the integrated radial 
momentum equation becomes 



dV r dV r GM C 2 <9E J 2 B z B rs HB Z dB, 



at 



8r 



E dr 



H d 



2ttE 

„2 D 2 



2ttE dr 



1 



+ m& dr W - ^ -6^Uj {B " ~ * 
When integrating the angular momentum equation (A3), we write 

d 



B 7 dBj, B r d , \ 1 d , \ B r d , N B^ . 
4tt & 4irrdr y <PJ 4ir dz V 1 07 4th- <9r V \itr dr r ' 



and use the same ansatz as above for B r (r,z) and B^{r,z). We then obtain 



dJ +v dJ 

dt ' dr 



rB.B 



(p 7 S 



H d 



{r 2 B ryS B 4 



rBrsB^s ( dH 
dr 



(A9) 



(A10) 



(All) 



27rE bnrY, dr x ''" ' v '" ' 37rE 
In integrating the vertical hydrostatic-balance equation, we take the pressure at the disk surfaces to vanish. The midplane pressure 



p c w {T I /2H)C 2 is then given by 



7T 2 GM C Y,H B 2 S + B^. HB r ,dB z 

^=2 GS + ^^ + — ^ 8^Ta7' 



(A12) 



where the first two terms on the right-hand side constitute our approximation to g z and represent, respectively, the local disk self- 
gravity and the tidal squeezing by the central point mass. 
The ion equation of motion is approximated by 

— V = -^(VxB)xB, (A13) 

where Vd = V, • — V„; we assume that the disk is sufficiently weakly ionized that the neutral velocity V„ is practically indistinguishable 
from the bulk velocity V. The radial component of this equation is 



4irp 

Tni 



dB r dB. 



Ba, d 



which yields, upon integration over z, 

Tni 



Vn r - 



2ttE 



, dB z 

B z I B rs — H-^- 



6ttE 



H 

— B rs 



dr 

dB 



r dr 



dr 



— Ba 



dB 



<p,S 



dr 



b: 



\ ( dH 

s > \~dr~ 



(A14) 



(A15) 
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(see eq. [A9]). Similarly, the azimuthal component of equation (A13), 

4-np 



-V, 



gives, when integrated over z, 



Vu = Vi + 



D.d 



dBj, B r d , s 

B ^ + Td~A rB ^ 



«h d 



V B r,sB<l>,s) - T-^ B r,sB^r — 



dH\ 
~dr~) 



(A16) 



(A17) 



6^rdr y -'<-v*j 3?rE - 
(see eq. [All]). 

The flux conservation equation is obtained from Faraday's law (dB/dt = — cV x E) and Ohm's law in the ambipolar diffusion- 
dominated regime (cE = —V,- x B) together with equation (A13). The result is 



?=Vx(VxB) + Vx(-^ 
dt {4-Kp 



[(V x B) x B] x B 



(A18) 



The z component of this equation is 

as* 
dt 



1 d 

2irr dr 
ld_ 
r dr 

-B r 



B, 



~dt) 

-rV r B z -\ 

dB^ 
' dz ' 



4np \ 
dB r 
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dr 
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which implies 
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(A20) 



Integrating this over z and rearranging, gives 

a* 

~dt ~ ~ 
1 

+ 



2irrV r B z + '■ 



— B rs —H- 



dB 7 



-B rs +HB Z 



— ln(rB^) 
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rBi 
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HB 7 
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(A21) 



In formulating the self-similar model in § 2.2 we use a pared-down version of the above equations. This is motivated by our finding 
that our basic results can be derived without including any O (H jr) terms. The only instance in which a term of this order plays a role 
in our solutions is in allowing us to refine the structure of the ambipolar-diffusion shock by employing the radial gradient of B z in the 
radial momentum equation. We therefore retain the combination [B r , s — H(dB z /dr)] (which is proportional to the vertically integrated 
azimuthal component of the current density) in equation (A9) but neglect the other two O(Hfr) terms in that equation, since they 
typically do not exceed the oc dB 7 jdr term. [The last two terms in equation (A9) can become much larger than the oc dB 7 jdr term in 
the limit r/t — > of the IMHD solution, when B r s /B z formally diverges. However, as we argue in § 3, this divergence is not expected 
to occur in a real disk. Furthermore, the neglect of these terms is justified even in our formal solution since H oc 1 /B 2 rs in this limit, 
so the radial derivative terms remain negligible in comparison with the magnetic tension (ex B z B r s ) term.] Since the same integrated 
current density component appears also in equations (A15) and (A21) for Vo. r and d^/dt, respectively, we keep the combination 
[B ns —H(dB z / dr)] in these equations too for self-consistency, but we again omit all the other O(Hfr) terms. For the same reason we 
keep the oc H(dB z /dr) term also in the vertical force balance equation (A12), which means that we end up retaining all the 0(H/r) 
terms in this "master equation" for H. 

We similarly neglect the 0(H/r) terms in the angular momentum equation (Al 1) and, correspondingly, in equation (A17) for V^. 
The latter equation can then be substituted into equation (3) to yield a purely algebraic expression for B^. 



B. SINGULAR LINES AND SHOCK JUMP CONDITIONS 

The singular lines of the flows considered in this paper can be obtained by rewriting the radial momentum equation (21) in the 
form 



da 



db 7 



(l-w z ) — = -b 7 h—- 
dx " dx 



(Bl) 



In the AD limit, in which the field is effectively decoupled from the bulk of the matter, we immediately infer that the singular line 
corresponds to the sonic speed: 

(x-u) 2 =l (AD limit) (B2) 

(see CCK for a discussion of the physical basis of this result). 

In the ideal-MHD limit, flux freezing implies b z j a = l/po (see eq. [33]), and therefore the singular line is given by 



(x - uf = 1 + ha/pl (IMHD limit) . 



(B3) 
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In the region where self-gravity dominates, h w 2/<r, and equation (B3) reduces to (x — u) 2 = 1 + 2//ip, which is the result given by 
equation (41) in CCK. 5 In dimensional form, equation (B3) is 

(r/t - V r ) 2 = C 2 +B?///2ttS = C 2 +B 2 /4tt P , (B4) 

where the last term on the right-hand side is equal to the Alfven speed in the disk. The IMHD singular line thus corresponds to the 
fast-magnetosonic speed, which agrees with (and generalizes) the statement made in CCK on the basis of the self-gravity-dominated 
limit of equation (B3). 

When the curve describing the flow crosses and falls below the singular line in the {x, u} plane, this generally signals the appearance 
of an unresolved shock discontinuity. In the present work, this is the situation invariably encountered in the vicinity of the centrifugal 
radius. To properly model this discontinuity, appropriate shock jump conditions need to be applied. The conservation of mass across 
the shock yields 

aw s = const, (B5) 

where w s = x s — u (with x s denoting the shock location). When the ions and neutrals are well coupled inside the shock, the magnetic- 
flux conservation condition takes the form 

b z w s = const . (B6) 

In the strong-coupling case, equations (B5) and (B6) imply fi = a/b z = const across the shock: in what follows we set \i = /i , 
which applies when the entire flow upstream of the shock obeys ideal MHD. Under the adopted isothermal approximation, the jump 
conditions are fully specified once we integrate also the differential terms in equation (Bl) across the discontinuity. Performing this 
integral requires making some assumptions about the behavior of a, b z , and h in the shock region. In the following subsections, 
we consider four cases in which algebraic jump conditions can be derived. These results have enabled us to explicitly model the 
centrifugal shock in all the core-collapse solutions presented in this paper. 

B. 1 . Isothermal Shock 

If one can neglect variations in the magnetic field across the discontinuity, as is the case when the centrifugal shock occurs deep 
inside the AD regime (see § 3.3.1), equations (Bl) and (B5) imply 

a{w] + 1) = const. (B7) 

Denoting the upstream and downstream sides of the shock by the subscripts 1 and 2, respectively, one can use equations (B7) and 
(B5) to express the nontrivial {w s \ ^ w S 2) solution in the form w s \w s2 = 1. This is the simplest realization of the "isothermal" jump 
conditions discussed by Shu & Li (1997). 

B.2. Generalized Isothermal Shock 

If the shock occurs in the self-gravitating regime, where h w 2/ a, then b z h oc b z j a is constant across the shock, and the integration 
of equation (Bl) leads to 

a(w 2 + &) = const, (B8) 

where 0=1 +2//1q. This situation is encountered in the fast-rotation solutions presented in § 3.2.2 and § 3.3.2. In this case the 
upstream and downstream radial velocities (in the shock frame) are related by w s iw S 2 = 0, which constitutes another realization 
of the "isothermal" jump conditions, remark. [Note that, as w S 2 — > w s i, one recovers the singular-line relation (eq. [B3]) for a 
self-gravitating IMHD flow.] 

B.3. Magnetically Squeezed Shock 

If magnetic squeezing due to the radial field dominates gravity in the vertical force-balance equation, then h w 2abj} (see eq. 
[42]), and equation (Bl) yields 

2er 3 

a(w 2 + l)-\ tpt= const. (B9) 



Taking into account the relation (B5), we obtain a quartic equation for er 2 (whose coefficients are combinations of the upstream flow 
variables). We can eliminate the trivial solution a 2 = &\ by dividing this equation by (02 — 01 ), which reduces it to 

a\ + aia 2 + (of + Q) a 2 - a x {u x -x s ) 2 Q = , (BIO) 

where Q = 3^lb 2 s /2. Equation (BIO) is a cubic of the form a\ +a2<7 2 +a\(T2 +«o = 0, whose discriminant is D = {p/'if + (q/2) 2 , 
with p = a\— a\j1 and q = «o — a^jT* + 2a\j21 . For this particular cubic, the discriminant is always positive, guaranteeing that 
there is exactly one real root, given by 

02 = -a 2 /3 +q+ + q- , (Bll) 

where q± = (—q/2± ^f~D) 1 ^. This solution provides a good approximation to the centrifugal shock in the fiducial IMHD flow 
(§3.2.1). 

5 Equation (41) in CCK is therefore only a special case of the more general expression, which is given by equation (B3). 
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B.4. Constant-Thickness Shock 

Another possibility is that h is constant across the shock, as would be the case if the central point-mass gravity controlled the 
squeezing in the vertical equilibrium equation (see eqs. [57] and [71]). In this case, integration of equation (Bl) results in 

cr(w 2 + l)+b 2 h/2 = const. (B12) 

In conjunction with the mass conservation equation (B5), this can be written as a cubic equation for a^. After factoring out the trivial 
solution, one obtains a quadratic equation, 

aj + (ai +2fio/h)a 2 -(m/x s ) 2 (2fi Q /h<T l )=0 , (B13) 

which has exactly one positive real root. In order for this result to be applicable to a centrifugal shock in a collapsing magnetized 
core, the shock must occur far enough from the center that ideal MHD is still a good approximation (so that the magnetic term in eq. 
[Bl] continues to play a role). However, in the solutions explored in this paper, the requirements that the central gravity dominate the 
vertical squeezing and that the IMHD regime be applicable are never satisfied at the location of the centrifugal shock: in the fiducial 
and strong-gravity IMHD cases, magnetic squeezing exceeds the tidal force in the central gravity-dominated region, whereas in the 
fast-rotation IMHD and AD cases, self-gravity dominates at the location of the shock. 

C. INCORPORATION OF A CENTRIFUGALLY DRIVEN WIND INTO THE DISK SOLUTION 

As we note in § 4, the asymptotic ambipolar-diffusion disk solution satisfies the launching condition B rs /B z > 1/V3 for a cold, 
centrifugally driven wind, and, in fact, implies a radial scaling of the magnetic field components (at any given instant of time) that 
is identical to that of the time-independent, radially self-similar wind solution of Blandford & Payne (1982, hereafter BP). Once a 
super- Alfvenic wind of this type is established, it constitutes the dominant mechanism of removing angular momentum from the disk 
(at least in the vertical direction). Furthermore, if the outflow is quasi steady on the accretion time scale, then the spatially self-similar 
wind solution can be used to determine the value of B^ (or, equivalently, of the parameter 8) in the asymptotic self-similar (in space 
and time) disk solution. An explicit wind solution also makes it possible to evaluate the reduction in the mass accretion rate onto 
the central star caused by the diversion of part of the inflowing mass into a disk outflow. Motivated by these considerations, we now 
show how the centrifugally driven wind solution can be incorporated into the asymptotic AD disk model. As demonstrated in § 3.3.1, 
the numerical solution of the disk equations typically converges rapidly behind the centrifugal shock to the asymptotic form given by 
equations (51)— (57). This indicates that, if the model parameters do not strongly differ from the fiducial values, then the results that 
we obtain for the asymptotic regime should be applicable to the bulk of the rotationally supported disk that forms around the central 
object. 

The mass loss in the BP model is constant for each decade of disk radius and can be related to the accretion rate M c at the inner 
edge of the disk by 

1 dM w = eM c 1 
2irr dr 2A - 3 2irr 2 ' ( ' 

where e ^ 1 is the fraction of the disk binding energy at its inner edge r m that is carried off by the wind, and where A is the total 
(kinetic plus magnetic) specific angular momentum in the wind, normalized by the Keplerian disk value rVx — [GM c r) x l 2 . The 
field-line constant A is one of the three parameters that define a solution in the cold, self-similar BP wind model, the other two being 
£q = B v /B z and k = 4iypV z Vk/B 2 (the normalized mass-to-magnetic flux ratio evaluated at the disk surface). As we noted above, for 
any given value of t in the asymptotic AD disk solution, B rs /B z is a spatial constant (=4/3; see eq. [56]) that exceeds the launching 
threshold 1 /\/3. It is thus natural to set £q = 4/3. One can then evaluate the parameter k by equating 2pV z with the expression 
given by equation (CI) (where the factor of 2 accounts for the two disk surfaces). The result (in terms of the nondimensional flow 
quantities introduced in § 2.3) is 

em c m c ' ,„„ N 

(C2) 



' (2X-3)b 2 z x 5 / 2 • 

Substituting for m c and b z from the asymptotic-solution equations (51) and (55), respectively, this becomes 

_ 2Se 
K ~2A^3 ' 

We now substitute for <5 in the asymptotic disk solution the corresponding expression in the BP wind model, 



(C3) 



5 = -^=k{\-\), (C4) 

which yields 



Although we started by evaluating k, we ended up with an expression that gives A in terms of e. The latter parameter was treated by 
BP as being distinct from A, but we are now able to relate them by combining the disk and wind solutions. By fixing the values of 
£q and A, one can derive the value of k from the solution of the self-similar wind equations (see Fig. 2 in BP), and thereby obtain S 
(using eq. [C4]). In practice, this "recipe" for incorporating the wind model does not alter our basic disk model since it effectively 
amounts to replacing one somewhat arbitrary parameter (<5) by another (e). However, by combining the disk and wind solutions, it is 
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possible to constrain the parameter values. In particular, by requiring that the flow be magnetically dominated in the vicinity of the 
disk (which implies n < 1) and that it become super- Alfvenic at a finite height, one can infer a lower bound on A. Based on Figure 2 
in BP, this lower limit is ~ 3 for £q = 4/3 (which, by eq. [C5], implies that e must be > 0.75), and as A increases from <~ 3 to <~ 30 
(corresponding to e increasing from ~ 0.75 to ~ 0.98), 5 decreases from ~ 2 to <~ 0.3. 

In the absence of a disk outflow, the central mass eigenvalue m c is given approximately by the product — xua evaluated at x c (see 
eqs. [26] and [51]). When the mass loss between the outer edge of the disk (r out w r c ) and its inner edge (r m w the protostellar 
radius) is taken into account in the mass conservation equation, the value of m c is reduced to 

(— xua) Xc , (C6) 

where A = r out /V; n (see eqs. [A5], [CI], and [C5]). To the extent that deuterium burning causes the radius of the accreting protostar 
to increase roughly linearly with its mass (e.g., Stahler 1988), A can be approximated as a constant. Estimating r c from our fiducial 
AD solution (see eq. [65]) and R* from the mass-radius relation given by Stahler (1988), we obtain In A w 8. Equation (C6) then 
implies that mass outflow from the disk surfaces could reduce the accretion rate onto the central object by a factor as large as <~ 3 
(using the minimum value of A, which corresponds to k = 1). More strongly magnetized outflows have lower m c reduction factors: 
for example, when k decreases to < 0.1 (corresponding to A > 10 and 5 < 1), this factor declines to < 1.5. Although the cumulative 
effect of the mass loss from the wind surfaces (embodied in the factor In A) can be significant, its impact at any given radius is 
comparatively small because of the large range of radii involved (r out /r; n «3x 10 3 in the fiducial case). The incorporation of a 
disk wind modifies the integrated self-similar disk equations only through logarithmic terms in x, so to lowest order the power-law 
scalings of the asymptotic disk solution remain unchanged. 

Explicit models of AD-dominated, wind-driving accretion disks were previously constructed by Wardle & Konigl (1993) and Li 
(1996). Two general types of solutions were obtained, depending on the value of the neutral-ion coupling parameter ?7wkl = 1 /t^^k- 
Strongly coupled disks have t/wkl > 1 at all heights and are characterized by the thermal pressure not being much larger than the 
magnetic pressure at z = 0, by B r starting to increase already near the midplane and generally exceeding \B$\ at the disk surface 
(which leads to a strong magnetic squeezing of the disk), and by the inflow speed typically being > C (Wardle & Konigl 1993). 
In contrast, in weakly coupled disks (characterized by rywKL < 1 over the bulk of the vertical column) the thermal pressure can be 
much larger than the magnetic pressure at z = 0, B r only starts to grow well above the midplane (at a height where the coupling 
parameter finally rises above 1) and usually does not exceed \B^,\ at the surface (which results in the magnetic squeezing remaining 
relatively unimportant), and the mass-averaged inflow speed is typically well below C (Li 1996). As we point out in § 4, 7/wkl is 
equal to 2/3(5 in the asymptotic AD disk solution. This result is compatible with the just-described behavior of B rs and |B^ jS |: since 
B r , s /B z = 4/3 in the asymptotic solution and 5 = \B,f, tS \/B z , the expression for ?7wkl implies B r ^/\B,f,^\ — 2?7wkl [consistent with eq. 
(4.5) in Wardle & Konigl (1993) and eq. (39) in Li (1996)]. The values of S derived from the BP wind solution indicate that our AD 
disks are likely to be weakly coupled. This conclusion is consistent with the low magnetic-to-thermal pressure ratio implied by the 
asymptotic disk solution (it is oc r 1 / 2 as r — > at any given instant of time; see eqs. [54]-[57]), which corresponds to the situation 
in weakly coupled disks. The association with the weakly coupled class of solutions is further indicated by the unimportance of 
magnetic squeezing in the asymptotic expression for the disk half-thickness (eq. [57]) and by the subsonic inflow speed inferred in 
the asymptotic regime (V r /C oc r 1 / 2 as r — > and t is held fixed; see eq. [53]). The inner regions of real protostellar disks (on radial 
scales ~ 0.1-10AU) are, in fact, expected to be so weakly ionized over most of their interiors that the requisite field-line bending 
for launching a centrifugal wind can only occur in comparatively thin surface layers (e.g., Wardle 1997). As was, however, already 
emphasized by Li (1996), the properties of such disks and of their associated outflows depend sensitively on the details of the spatial 
profile of ?7wkl- Further progress in developing these models will thus require a good understanding of the ionization structure in the 
inner regions of protostellar disks. 
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